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ABSTRACT 


This  work  addresses  the  design  and  analysis  of  the  discrete  longitudinal 
autopilot  for  application  to  the  bank-to-turn  (BTT)  missiles. 

In  the  development  the  classical  design  and  analysis  were  reviewed 
using  the  continuous  uncoupled  pitch  channel  autopilot  for  circular  airframe. 
Then,app1ying  analog-to-digital  conversion,  the  corresponding  discrete 
autopilot  was  designed.  The  performance  of  both  continuous  and  discrete 
open  loop  systems  was  analyzed  according  to  the  desired  requirements. 

In  the  following  section,  utilizing  modern  control  design  techniques  for 
the  discrete  pitch  channel  autopilot,  the  control  law  was  designed, 
assuming  availability  of  all  the  states  and  taking  the  desired  poles  to  be  at 
the  same  place  as  in  the  classical  desiyn.  Therefore  the  discrete  state- 
feedback  autopilot  was  obtained  and  analyzed.  The  next  step  was  the  design 
of  an  estimator  which  estimates  the  entire  state  vector,  given 
measurements  of  the  portion  of  the  state.  The  performance  of  the  estimator 
was  analyzed  too. 

Finally  coupling  the  discrete  pitch  and  roll  channel  autopilots  the 
overall  performance  of  the  discrete  system  was  obtained  and  analyzed.  The 
resulting  design  was  also  proved  to  be  robust.  


TABLE  OF  CONTENTS 

I.  INTRODUCTION 19 

II.  CLASSICAL  LINEAR  DESIGN  AND  ANALYSIS  OF 
LONGITUDINAL  UNCOUPLED  CHANNEL  AUTOPILOT  FOR 
CIRCULAR  AIRFRAME 23 

A.  GENERAL 23 

B.  AIRFRAME  CONFIGURATION 24 

C.  UNCOUPLED  LINEAR  PITCH  CHANNEL  AUTOPILOT 26 

1 .  Aerodynamic  Model 26 

2.  Pitch  Control  Law 28 

3.  Transfer  Functions  of  Aerodynamic  Model 33 

4.  Design  Approach  and  Analysis  of  the 
Continuous  Open  Loop  System 34 

5.  Design  Approach  and  Analysis  of  the 

Discrete  Open  Loop  System 56 

6.  Comparison  of  Continuous  and  Discrete 

Systems 65 

III.  MODERN  CONTROL  DESIGN  AND  ANALYSIS  OF 
LONGITUDINAL  UNCOUPLED  CHANNEL    AUTOPILOT    66 

A.  GENERAL 66 

B.  STATE  FEEDBACK  DESIGNED  AUTOPILOT 68 

1.  Design  Approach  and  Analysis  of 

Control  Law 68 

2.  Performance  of  System 70 


3.     Simplified  State  Feedback  Designed 

Autopilot 78 

C.      ESTIMATOR  DESIGNED  AUTOPILOT 79 

1.  Design  Approach  and  Analysis  of 

Discrete  Closed  Loop  Estimator 79 

2.  Performance  of  System 88 

3.  Simplified  Estimator  Designed  Autopilot 97 

IV.  DISCRETE      COUPLED      PITCH      AND      ROLL      CHANNEL 
AUTOPILOTS    105 

A.  GENERAL 105 

B.  DESIGN  APPROACH  AND  ANALYSIS  OF  THE 

DISCRETE  STATE    FEEDBACK  AUTOPILOT 106 

C.  PERFORMANCE  OF  SYSTEM 106 

D.  ROBUSTNESS  OF  SYSTEM 1 14 

V.  CONCLUSIONS  AND-  RECOMMENDATIONS  122 

APPENDIX    A:      GEOMETRIC  AND  MASS  PROPERTIES  OF 

MISSILE  CONFIGURATION 123 

APPENDIX  B:  ANALOG  TO  DIGITAL  CONVERSION  OF 

UNCOUPLED  PITCH  CHANNEL  AUTOPILOT 124 

APPENDIX   C:      PROGRAM  LOGIC  FOR  COMPUTING  CONTROL 

LAW  K  VIA  ACKERMANN'S  FORMULA 125 

APPENDIX   D:      COMPUTER  PROGRAM  FOR  APPLICATION  OF 

ACKERMANN'S  FORMULA 127 

APPENDIX   E:      ANALOG  TO  DIGITAL  CONVERSION  OF 

COUPLED  PITCH  AND  ROLL  CHANNEL 


AUTOPILOTS 129 

APPENDIX   F:      ROBUSTNESS  DATA  FILE  FOR  THE  DISCRETE 

STATE  FEEDBACK  DESIGNED  AUTOPILOT 131 

APPENDIX   G:      ROBUSTNESS  DATA  FILE  FOR  THE  DISCRETE 

ESTIMATOR  DESIGNED  AUTOPILOT 134 

LIST  OF  REFERENCES    1 37 

INITIAL  DISTRIBUTION  LIST 1 38 


LIST  OF  TABLES 

I.  Linearized  Aerodynamic  Derivatives  (M=3.95) 35 

II.  Stability  Margins  of  Uncoupled  Pitch  Channel 51 


8 


LIST  OF  FIGURES 

2.1  Bank-to-Turn  Autopilot .25 

2.2  Model  of  Circular  Configuration- 1  /6  Scale 27 

2.3  Aerodynamic  Sign  Convention 29 

2.4  Linear  Pitch  Channel  Dynamic  Model 30 

2.5  Pitch  Control  Law 31 

2.6  Uncoupled  Pitch  Channel  Autopilot 32 

2.7  Pitch  Normal  Acceleration  vs  Time;  Uncoupled  Pitch 
Channel;  Continuous  Open  Loop  System 42 

2.8  Pitch  Angular  Rate  vs  Time;  Uncoupled  Pitch  Channel- 
Continuous  Open  Loop  System 43 

2.9  Pitch  Tall  Incidence  vs  Time;  Ucoupled  Pitch  Channel; 
Continuous  Open  Loop  System 44 

2.10  Pitch  Normal  Acceleration;  Uncoupled  Pitch  Channel; 
Continuous  Open  Loop  System;  Gain  vs  Frequency 45 

2.1 1  Pitch  Normal  Acceleration;  Uncoupled  Pitch  Channel; 
Continuous  Open  Loop  System;  Phase  vs  Frequency 46 

2. 1 2  Pitch  Angular  Rate;  Uncoupled  Pitch  Channel; 

Continuous  Open  Loop  System;  Gain  vs  Frequency 47 

2. 1 3  Pitch  Angular  Rate;  Uncoupled  Pitch  Channel; 

Continuous  Open  Loop  System;  Phase  vs  Frequency 48 

2. 1 4  Pitch  Tail  Incidence;  Uncoupled  Pitch  Channel; 

Continuous  Open  Loop  System;  Gain  vs  Frequency 49 

2.15  Pitch  Tail  Incidence;  Uncoupled  Pitch  Channel; 


Continuous  Open  Loop  System;  Phase  vs  Frequency 50 

2.16  Pitch  Normal  Acceleration;  Uncoupled  Pitch  Channel; 
Continuous  Open  Loop  System;  Pole-Zero  Map 52 

2. 1 7  Pitch  Angular  Rate;  Uncoupled  Pitch  Channel; 

Continuous  Open  Loop  System;  Pole-Zero  Map 53 

2.18  Pitch  Tail  Incidence;  Uncoupled  Pitch  Channel; 

Continuous  Open  Loop  System;  Pole-Zero  Map 54 

2. 1 9  Pitch  Normal  Acceleration  vs  Time;  Uncoupled  Pitch 
Channel;  Discrete  Open  Loop  System 57 

2.20  Pitch  Angular  Rate  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Open  Loop  System 59 

2.2 1  Pitch  Tail  Incidence  vs  Time;  Ucoupled  Pitch  Channel; 
Discrete  Open  Loop  System 60 

2.22  Pitch  Normal  Acceleration;  Uncoupled  Pitch  Channel; 
Discrete  Open  Loop  System;  Pole-Zero  Map 61 

2.23  Pitch  Angular  Rate;  Uncoupled  Pitch  Channel; 

Discrete  Open  Loop  System;  Pole-Zero  Map 62 

2.24  Pitch  Tail  Incidence;  Uncoupled  Pitch  Channel; 

Discrete  Open  Loop  System;  Pole-Zero  Map 63 

3.1  Pitch  Normal  Acceleration  vs  Time;  Uncoupled  Pitch 
Channel;  Discrete  State-Feedback  Autopilot 72 

3.2  Pitch  Angular  Rate  vs  Time;  Uncoupled  Pitch  Channel- 
Discrete  State-Feedback  Autopilot 73 

3.3  Pitch  Tail  Incidence  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  State-Feedback  Autopilot 74 


10 


3.4  Pitch  Normal  Acceleration  ;  Uncoupled  Pitch  Channel; 
Discrete  State-Feedback  Autopilot;  Pole-Zero  Map 75 

3.5  Pitch  Angular  Rate;  Uncoupled  Pitch  Channel;  Discrete 
State-Feedback  Autopilot;  Pole-Zero  Nap 76 

3.6  Pitch  Tall  Incidence;  Uncoupled  Pitch  Channel;  Discrete 
State-Feedback  Autopilot;  Pole-Zero  Map 77 

3.7  Pitch  Normal  Acceleration  vs  Time;  Uncoupled  Pitch 
Channel;  Discrete  Simplified  State-Feedback  Autopilot 80 

3.8  Pitch  Angular  Rate  vs  Time;  Uncoupled  Pitch  Channel- 
Discrete  Simplified  State-Feedback  Autopilot 81 

3.9  Pitch  Tall  Incidence  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Simplified  State-Feedback  Autopilot 82 

3. 1 0  Pitch  Normal  Acceleration  ;  Uncoupled  Pitch  Channel; 
Discrete  Simplified  State-Feedback  Autopilot;  Pole- 
Zero  Map 83 

3. 1 1  Pitch  Angular  Rate;  Uncoupled  Pitch  Channel;  Discrete 
Simplified  State-Feedback  Autopilot;  Pole-Zero  Map 84 

3. 1 2  Pitch  Tail  Incidence;  Uncoupled  Pitch  Channel;  Discrete 
Simplified  State-Feedback  Autopilot;  Pole-Zero  Map 85 

3. 1 3  Open  Loop  Estimator 87 

3. 1 4  Closed  Loop  Estimator 87 

3. 1 5  Pitch  Normal  Acceleration  vs  Time;  Uncoupled  Pitch 
Channel;  Discrete  Estimator  Autopilot 90 

3. 1 6  Pitch  Angular  Rate  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Estimator  Autopilot 91 


11 


3. 1 7  Pitch  Tail  Incidence  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Estimator  Autopilot 92 

3.18  Pitch  Normal  Acceleration  ;  Uncoupled  Pitch  Channel; 
Discrete  State  Feedback  and  Estimator  Autopilots; 
Pole-Zero  Map 93 

3. 1 9  Pitch  Angular  Rate;  Uncoupled  Pitch  Channel;  Discrete 
State  Feedback  and  Estimator  Autopilots;  Pole-Zero 

Map 94 

3.20  Pitch  Tail  Incidence;  Uncoupled  Pitch  Channel;  Discrete 
State  Feedback  and  Estimator  Autopilots;  Pole-Zero 

Map 95 

3.2 1  Pitch  Normal  Acceleration  vs  Time;  Uncoupled  Pitch 
Channel;  Discrete  Simplified  Estimator  Autopilot 99 

3.22  Pitch  Angular  Rate  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Simplified  Estimator  Autopilot 100 

3.23  Pitch  Tall  Incidence  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Simplified  Estimator  Autopilot 101 

3.24  Pitch  Normal  Acceleration  ;  Uncoupled  Pitch  Channel; 
Discrete  Simplified  Estimator  Autopilot;  Pole-Zero 

Map 102 

3.25  Pitch  Angular  Rate;  Uncoupled  Pitch  Channel;  Discrete 
Simplified  Estimator  Autopilot;  Pole-Zero  Map. 103 

3.26  Pitch  Tall  Incidence;  Uncoupled  Pitch  Channel;  Discrete 
Simplified  Estimator  Autopilot;  Pole-Zero  Map 104 

4. 1        Pitch  Normal  Acceleration  vs  Time;  Coupled  Pitch  and 


12 


Roll  Channels;  Discrete  State-Feedback  Autopilot 108 

4.2  Pitch  Angular  Rate  vs  Time;  Coupled  Pitch  and  Roll 
Channels;  Discrete  State-Feedback  Autopilot 109 

4.3  Pitch  Tall  Incidence  vs  Time;  Coupled  Pitch  and  Roll 
Channels;D1screte  State-Feedback  Autopilot 1 10 

4.4  Pitch  Normal  Acceleration  ;  Coupled  Pitch  and  Roll 
Channels;  Discrete  State-Feedback  Autopilot;  Pole-Zero 

Map Ill 

4.5  Pitch  Angular  Rate;  Coupled  Pitch  and  Roll  Channels; 
Discrete  State-Feedback  Autopilot;  Pole-Zero  Map 1 12 

4.6  Pitch  Tail  Incidence;  Coupled  Pitch  and  Roll  Channels; 
Discrete  State-Feedback  Autopilot;  Pole-Zero  Map 1 13 

4.7  SVADMO  vs  Frequency;  Coupled  Pitch  and  Roll  Channels; 
Discrete  State-Feedback  Autopilot 1 16 

4.8  SVADMO  vs  Frequency;  Coupled  Pitch  and  Roll  Channels; 
Discrete  Estimator  Autopilot 1 17 

4.9  MIN.ADD.IN.SV  vs  Frequency;  Coupled  Pitch  and  Roll 
Channels;  Discrete  State-Feedback  Autopilot 1 18 

4.10  MIN.ADD.IN.SV  vs  Frequency;  Coupled  Pitch  and  Roll 
Channels;  Discrete  Estimator  Autopilot 1 19 


13 


TABLES  OF  SYMBOLS  AND  ABBREVIATIONS 

BTT  Bank-to-Turn 

CBTT  Coordinated  Bank-to-Turn,  minimum  sideslip,  positive  a, 

<P«il80deg. 

Cy  rolling  moment  coefficient 

Cm  pitching  moment  coefficient 

C™  slope  of  curve  of  pitching  moment  coefficient  Cmvsa 

Cmsp  change  in  C™  per  degree  pitch  control  incidence,^ 

CN  normal  force  coefficient 

C^  slope  of  curve  of  normal  force  coefficient  CNvsa 

Cflfip  change  in  C^  per  degree  pitch  control  incidence,  6p 

CY  side  force  coefficient 

d  reference  length  for  coefficients  (=  2  ft) 

Iyy  moment  of  inertia  about  &  axis 

Izz  moment  of  inertia  about  z^  axis 

Ixx  moment  of  inertia  about  x^  axis 

KA  autopilot  pitch  acceleration  error  gain 

KYp  CBTT  autopilot  coordination  branch  gain 

p  roll  rate  about    x^ 

p  roll  acceleration  about    x^ 

POC  preferred  orientation  control 

q~  dynamic  pressure 

q  pitch  rate  about  y^ 

q  pitch  angular  accelaration  about  j& 


14 


Q,  constant  or  equilibrium  pitch  angular  rate 

r  yaw  angular  rate  about  z^ 

r  yaw  angular  acceleration  about  z^ 

S  reference  area  for  coefficients  (=fl  ft2) 

STT  Skid-to-Turn,  roll  attitude  stabilized 

u  velocity  component  in   x^  direction 

v  velocity  component  in  y^  direction,  assumed  to  constant 

V  constant  missile  flight  path  velocity 

V  missile  velocity  vector 

w  velocity  component  in    z^  direction 

x^  body-fixed  roll  axis,  along  axis  of  symmetry,  positive  forward 

y^  body-fixed  pitch  axis,  positive  starboard 

z^  body-fixed  yaw  axis,  forms  right  handed  orthogonal  system 

with  x^  and   & 

Ty  vehicle  axis  in  downward  direction  along  local  gravity 

vector ,  approximated  as  inertial  axis 

r±  achieved  normal  acceleration  in    z^  direction 

n20  commanded  normal  acceleration  in    TB  direction 

ny  achieved  normal   acceleration  in  y^  direction 

n?  achieved  normal   acceleration  in  Ty  direction 

ny  achieved  normal   acceleration  in  y^  direction 

lie  normal  acceleration  command  from  guidance  computer  in 

Zy  direction  plus  anti-gravity  bias  command 

nZc  normal  acceleration  guidance  command  in  z^  direction 

nYc  normal  acceleration  guidance  command  in  y^  direction 


15 


<pc  roll  attitude  command  from  guidance  computer,  zero  degrees 

in  -Ty  direction  and  90  degrees  in  y^  direction 
(p  roll  attitude,  zero  degrees  in  -%  direction  and  90  degrees 

in  y^  direction 
<pe  roll  attitude  error,  <pc-<p 

9  Elevation  Euler  Angle,  second  rotation,    J(q  coscp-rsin<p)dt 

Y  Azimuth  Euler  Angle,  first  rotation  about  y^  J(q  cos<p+rsin(p)dt 

5p  pitch  control  incidence  (positive  tail  incidence  produces 

negative  pitching  moment) 
5^  commanded  pitch  control  incidence,  6p 

$y  yaw  control  incidence  (positive  tail  incidence  produces 

negative  yawing  moment) 
5Yc  commanded  yaw  control  incidence,  Sy 

Sr  roll  control  Incidence  (positive  tail  incidence  produces 

positive  rolling  moment) 
Src  commanded  roll  control  incidence  ,  6r 

a*  constant  or  equilibrium  angle-of-attack 

a  angle-of-attack 

a  angle-of-attack  rate 

£  angle   of   sideslip 

£  sideslip  angular  rate 

OPTSYS  Optimal  Systems  Control  Program 

ORACLS  Optimal  Regulator  Algorithms  for  the  Control  of  Linear 

Systems 
X  plant  state  =  Ns  x  1  or  n  x  1 


16 


Ns 

number  of  states 

U 

control  =  1%  x  1  or  m  x  1 

No 

number  of  controls 

F 

continuous  plant  system  matrix 

A 

discrete  plant  system  matrix 

G 

continuous  plant  control  input  matrix 

B 

dicrete  plant  control  input  matrix 

y 

output  measurements 

H 

plant  output  matrix 

J 

plant  direct  transmission  matrix 

K 

control  gain  (law) 

L 

estimator  gain 

T 

sampling  period 

17 


ACKNOWLEDGEMENT 

I  would  like  to  express  my  gratitude  to  Professor  DJ.  Collins,  for  his 
valuable  guidance  and  assistance  that  helped  me  finish  this  work. 

I  would  also  like  to  dedicate  this  thesis  to  my  wife  Stavroula,  and  my 
parents  Christos  and  loanna  for  their  love  and  understanding. 


18 


I.      INTRODUCTION 


Missile  systems  designed  to  meet  the  threats  of  the  future  will  require 
capabilities  far  beyond  of  those  currently  available.  Future  requirements 
dictate  the  need  for  dramatic  increases  in  Short  Range  Air-to-Air  Missile 
(SRAAM)  maneuverability  and  guidance  accuracy.  Many  current  missile 
designs,  which  are  characterized  by  cruciform  airframes  and  skid-to-turn 
(STT)  steering,  cannot  be  extended  to  meet  these  requirements  due  to  low 
aerodynamic  efficiency.  Bank-to-turn  steering  provides  the  capability  to 
design  asymmetric  missile  airframes  which  optimize  aerodynamic 
performance  in  one  plane  and  thus  obtain  the  SRAAM  maneuverability  and 
accuracy  requirements.  Furthermore  the  ramjet  engine  chin  inlet 
configuration  is  believed  to  provide  a  greater  range  capability  than  other 
inlet  configurations,  and  BBT  control  is  required  to  satisfy  the  sideslip 
constraints  imposed  by  chin  inlets. 

In  addition,  to  its  potential  advantages,  BBT  steering  introduces  some 
technical  concerns  which  must  be  carefully  evaluated.  For  example,  the 
methodology  for  designing  a  bank-to-turn  autopilot  is  not  well  developed. 
Such  a  design  must  take  into  account  the  aerodynamic  and  kinematic 
coupling  terms  as  well  as  allow  for  operation  a  low  signal  levels  (i.e,  small 
angle  of  attack)  when  the  preferred  roll  orientation  is  poorly  defined.  In 
addition,  the  coupling  of  body  motion  into  the  guidance  signals  is  another 
major  concern  for  BBT  systems.  Recent  analyses  have  indicated  that  some 
skid-to-turn  (STT)  systems  can  tolerate  a  limited  amount  of  radome- 
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Induced  instabilities  without  severe  performance  degradation  [Ref.U.  In  a 
BBT  system  a  coupling  loop  is  closed  through  roll  rate  as  well  as  pitch  and 
yaw  rates  [Ref.2l.  It  is  not  known  whether  BBT  systems  can  tolerate 
coupling  induced  instabilities.  Another  concern  is  the  interaction  of  BBT 
control  with  missile  functions  such  as  detection  (seeker),  guidance  signal 
processing  ^control  surface  effectiveness,  etc. 

All  the  above  concerns  must  be  investigated  before  BBT  steering  can  be 
considered  a  viable  method  to  control  high  performance  tactical  missiles. 

Many  missile  programs  were  initiated  during  the  past  decade  to 
improve  the  capability  of  steering  a  tactical  missile  via  BTT  control.  The 
results  have  greatly  advanced  the  understanding  of  the  various  missile 
subsystems. 

In  the  autopilot  area,  specifically,  many  types  of  autopilots  have  been 
found  which  force  the  missile  to  roll  or  bank  so  that  the  steering  maneuver 
occurs  with  the  missile  airframe  oriented  in  a  specific  or  preferred 
direction  with  respect  to  the  incoming  airsteam.  This  entire  class  of 
autopilots  referred  to  as  Preferred  Orientation  Control  (POC)  autopilots. 

The  main  criterion  for  the  selection  of  one  particular  type  of  autopilot 
is  based  upon  the  guidance,  airframe,  and  propulsion  system  requirements. 
Generally  missiles  with  either  one  or  two  planes  of  symmetry  use  a  POC 
autopilot  which  forces  the  missile  to  bank  in  order  to  turn  (BTT)  like  an 
aircraft,  and  if  this  motion  is  coordinated  it  is  usually  referred  to  as  CBTT. 

In  order  to  take  full  advantage  of  CBTT  control,  planar  airframes  have 
been  designed  to  increase  lifting  capability  in  one  direction  without  the 
weight  and  drag  penalty  associated  with  orthogonal  lifting  surfaces  [Ref.lJ. 


20 


These  airframe  have  aerodynamic  properties  which  have  the  potential  to 
enhance  CBTT  control. 

The  present  work  addresses  the  design  and  analysis  of  the  discrete 
pitch  channel  autopilot  for  application  to  the  bank-to-turn  (BTT)  missiles, 
with  circular  airframe  configuration.  The  circular  airframe  is  chosen 
because  it  is  stable  for  all  angles  of  attack. 

In  the  development  the  classical  design  and  analysis  of  the  uncoupled 
pitch  channel  autopilot  as  were  developed  In  reference  [Ref.3],  are  reviewed. 
Using  state-space  representation  a  tenth  order  system  is  formulated  and 
the  Step  Response,  Bode  and  Pole-Zero  Nap  plots  are  obtained  and  analyzed 
according  to  the  desired  requirements.  Then  the  continuous  system  is 
converted  into  the  corresponding  discrete  system  using  Analog  to  Digital 
conversion  from  ORACLS,  and  analyzed  according  to  the  Step-Response,and 
Pole-Zero  Map  plots  obtained.  A  comparison  of  the  plots  of  both  the 
continuous  and  discrete  systems  shows  no  significant  differences. 

Next,  utilizing  modern  control  design  techniques  and  assuming  that  we 
have  all  the  states  at  our  disposal  for  feedback  purposes,  the  control  law  is 
designed  using  the  Ackermann's  formula  (Ref.  6]  and  taking  the  desired  poles 
to  be  at  the  same  place  as  in  the  classical  case.  Using  the  control  law  gains 
the  discrete  state-feedback  autopilot  is  designed  and  the  Step  Response, and 
Pole-Zero  Map  plots  are  obtained  and  analyzed.  Then,  since  we  usually  know 
only  a  portion  of  the  state  variables,  an  estimator  is  introduced  in  order  to 
estimate  the  entire  state  vector  given  measurements  of  that  portion.  The 
poles  of  the  estimator  are  moved  to  the  left  of  the  original  poles  of  the 
system  in  order  to  get  faster  response.  Plots  as  above  are  obtained  and  an 
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analysis  is  made.  Next  the  designed  discrete  pitch  channel  autopilot  is 
coupled  with  the  discrete  roll  channel  autopilot  and  the  whole  discrete 
system  is  analysed  and  proved  to  be  robust.  Finally  conclusions  and 
recommendations  for  future  study  are  stated. 

The  analysis  in  all  the  above  cases  was  performed  using  the  existing  at 
NPS  Optimal  Systems  Control  Program  (OPTSYS)  for  the  continuous  system, 
Optimal  Regulator  Algorithms  for  the  Control  of  Linear  Systems  (ORACLS) 
for  the  discrete  system,  and  Pole  Placement  and  Robustness  Design  Program 
(POPLAR)  for  the  robustness  of  the  discrete  system. 
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II.    CLASSICAL  LINEAR  DESIGN  AND  ANALYSIS  OF  LONGITUDINAL  UNCOUPLED 
CHANNEL  AUTOPILOT  FOR  CIRCULAR  AIRFRAME 

A.  GENERAL 

The  Initial  phase  In  the  design  of  the  CBTT  autopilot  Involved  the  design 
of  Individual,  uncoupled  channels,  pitch,  yaw,  and  roll,  with  prescribed 
relationships  between  speeds  of  response  which  would  meet  the  CBTT 
requirements  when  coupled. 

The  design  technique  for  the  uncoupled  pitch  channel  was  classical,  using 
a  combination  of  Frequency  Response  and  Root  Locus  techniques  [Ref.  2],  to 
achieve  practical  bandwidths  and  In  turn  provide  the  range  of  required 
missile  body  angular  rates  and  control  motions.  In  addition,  the  resulting 
design  is  found  to  be  robust. 

The  whole  design  was  done  for  a  flight  condition  at  M=3.95  and 
Alt1tude=60  KFT. 

The  requirements  for  the  classical  design  technique,  as  stated  in 
[Ref.  21  are: 

1.  High  Freouencu  Attenuation  in  Actuator  Command  Branch 

>  15  dB  at  100  rad/sec  and  zero  angle-of-attack.  This  requirement 
limits  autopilot  speeds  of  response. 

2.  BfllflUffl  SiflMUlUl 

Gain  margins  >  6  dB,  Phase  margins  >  30  degrees  with  a  goal  of 
12  dB  and  50  degrees. 
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3.  Acceleration  Time  Response 

a.  63  percent  time  constant  of  0.5  seconds  for  a  step  command  of 
acceleration  at  the  flight  condition  of  interest  (M=3.95,  Altitude=60  KFT) 
and  small  angle  of  attack.  This  response  is  representative  of  a  tactical 
missile  of  this  size. 

b.  Overshoot  <  10  percent. 

c.  Zero  steady  state  error  in  acceleration  to  reduce  variations  of 
guidance  navigation  gain. 

In  this  chapter  the  airframe  configuration,  the  aerodynamic  model  and 
control  law  for  the  pitch  channel  are  reviewed.  The  transfer  functions  for 
circular  airframe  which  is  stable  for  all  angles  of  attack  [Ref.  2],  are 
derived.  An  analysis  for  the  uncoupled  pitch  channel,  as  regards  the 
acceleration  response,  the  body  angular  rate  and  control  surface  deflection, 
In  terms  of  time  resposnes  and  pole-zero  map  plots  for  both  the  continuous 
and  discrete  system  and  frequency  responses  for  the  continuous  system,  is 
made,  using  the  existing  at  NPS  Optimal  Systems  Control  Program  (OPTSYS) 
for  the  continuous  system,  and  Optimal  Regulator  Algorithms  for  the  Control 
of  Linear  Systems  (ORACLS)  for  the  discrete  system.  Finally  a  comparison  of 
the  continuous  and  discrete  system  by  analyzing  the  above  results  is  made. 

A  general  block  diagram  of  the  BTT  autopilot  is  shown  in  Figure  2.1. 

B.  AIRFRAME  CONFIGURATION 

An  investigation  has  been  conducted  by  NASA  [Ref.4]  to  compare  the 
experimental  aerodynamic  characteristics  of  a  low-drag  missile  concept 
with  a  body  of  circular  cross  section  to  one  with  a  body  of  3:1  elliptical 
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cross  section,  the  bodies  having  identical  cross  section  area  distributions. 
The  concepts  were  of  mono-wing  design  with  constant  wing  span.  Tail 
surfaces  were  located  flush  at  the  body  base  with  ±30  degrees  dihedral. 
Wind-tunnel  tests  were  performed  at  Mach  numbers  from  0.5  to  4.63  and  at 
angles  of  attack  from  about  -5°  to  +28  °. 

The  comparison  shows  no  significant  subsonic  normal  force  differences 
at  low  angles  of  attack;  however,  at  supersonic  speeds,  the  elliptical 
concept  increasingly  provides  greater  normal  force  up  to  Mach  2.5  to  3.0, 
beyond  which  an  incremental  increase  of  about  25  percent  holds  through  the 
angle  of  attack  range.  More  pronounced  nonlinearities  in  pitching  moment 
occur  at  subsonic  speed  for  the  elliptical  concept,  as  well  as  less 
longitudinal  stability  at  all  test  Mach  numbers.  However,  levels  of 
directional  and  lateral  stability  are  increased,  especially  at  the  higher 
angles  of  attack. 

This  work  uses  the  circular  airframe  configuration  which  is  shown  in 
Figure  2.2  and  it  is  taken  from  [Ref.  41.  The  circular  cross  sectional  body  has 
a  closure  ratio  At^/A,^  of  0.69  and  the  A™*  occurs  at  68  percent  body 

length. 

C.  UNCOUPLED  LINEAR  PITCH  CHANNEL  AUTOPILOT 
1.  Aerodynamic  Model 

The  following  addresses  the  linear  design  and   analysis   of   the 
continuous  uncoupled  pitch  channel  autopilot  for  the  circular  airframe. 
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A  linearized  aerodynamic  model  was  developed  for  stability  studies. 
The  method  used  is  an  extension  of  the  linearization  technique  used  for 
skid-to-turn  (STT)  aerodynamic  models. 

The  following  three  assumptions  were  made: 

a.  Plane  XB-FB  of  Figure  2.3  is  the  maneuver  plane. 

b.  Missile  is  trimmed  in  pitch  (i.e,  My=0,  at  fixed  values  of  a,  q,  and 

v- 

c.  Missile  roll  rate  is  constant. 

The  resulting  model  is  shown  in  block  diagram  form  in  Figure  2.4. 
Aerodynamic  stability  derivatives  for  M=3.95  are  provided  in  Table  I. 
2.  Pitch  Control  Law 

The  pitch  control  law  for  the  circular  airframe  is  shown  in  Figure  2.5 
as  given  in  [Ref.  2].  Lag-leads  were  used  to  prevent  guidance  noise 
saturation  problems.  The  rate  error  compensation  determines  the  high 
frequency  attenuation  and  was  used  to  minimize  the  effect  of  aerodynamic 
variations  on  acceleration  time  response.  The  acceleration  error 
compensation  determines  the  acceleration  time  response.  An  integrator  was 
used  in  the  acceleration  error  branch  to  satisfy  the  guidance  requirement  of 
zero  steady-state  error. 

A  normal  acceleration  command  (n^)  equal  to  1  gee  is  applied  to  the 

pitch  control  law  which  uses  measurements  of  missile  body  pitch  angular 
rate  (q)  and  pitch  normal  acceleration  (nz)  to  determine  the  required 

actuator  command  (6pc).  The  actuator  is  modeled  as  a  first  order  lag  at 
188.4  rad/sec  (30  HZ)  and  is  shown  in  Figure  2.6. 
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3.  Transfer  Functions  of  Aerodunamtc  Model 

Uncoupled  pitch  aerodynamic  transfer  functions  are 

Es 
+  1 


q  k(AE-BC)  k(AE-BC) 


deg/sec/deg  ( II.C.3-1) 


5p  -C  s2       Aks 

—  + +  i 

-C         -C 
Bs2 


+  1 


i\  AE-BC  AE-BC 


6p  -C  s2       Aks 

—  ♦ ♦  i 

-C         -C 


g's/deg         (II.C.3-2) 


q  S  q  S  (57.3)    q  S  d 

W  W  Iyy 

(57.3)  q"Sd  1845 

E  =  Cms  p  t       k  = 


Where: 

V=Ma=3.95  x  968.47=  3825.46  ft/sec 

S=  71  f  t2 

d=2  ft 

q=  1650      lb/ft2 

1^=  804      slug-ft2 

W=  2525      lbs 

Cm.  =  -006 

Cms  p=  -0-08 

CNa  =0.15 

CN5  p  =  0.04 
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As  referred  to  in  Appendix  A  and  Table  I. 

Then 

A  =0.30793 
B  =0.08211 
C  =  -44.33 1 62 
E  = -59.10883 
k  =  0.48229 

And  substituting  the  aforementioned  values  into  equations  (II.C.3-1) 
and  (II.C.3-2)  they  become: 


-1.3361  s- 0.159 


(1/sec)       (II.C.3-3) 


8p  22.545  x  1 0-3  S2  +  3.3634  x  1 0"*  s  +  1 


n2  1 06.074  x  1 0"3  s2  -  1 8.82 


5p  22.545  x  1 0-3  s*  +  3.3634  x  1 0"*  s  +  1 


(g's/rad)    (II.C.3-4) 


Equations  (II.C.3-3)  and  (II.C.3-4)  are  the  pitch  aerodynamic  transfer 
functions  for  the  pitch  angular  rate  (q)  about  the  YB  and  the  achieved 
maneuver  acceleration  r^  in  the  ZB  direction. 

4.  Design  Approach  and  Analysis  of  the  Continuous  Ooen  Loop  Sustem 
A  fixed  flight  condition  at  60  KFT  altitude  and  M=3.95  was  selected 
for  these  preliminary  performance  studies  of  circular  airframe.  This  flight 
condition  provides  a  sufficiently  low  dynamic  pressure,  so  that  missile 
maneuvers  will  result  in  large  enough  angles-of-attack,  to  exercise  side 
slip  control.  The  above  fixed  condition  was  selected,  in  order  to  reduce  the 
complexities  of  CBTT  autopilots. 
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a=0° 

Cy# 

-0.065 

^n£ 

-0.025 

C|0 

0 

CYSY 

0.021 

CnSY 

-0.050 

ClSY 

0 

CY8R 

0 

^n5R 

0 

C15R 

0.031 

Cn* 

0.15 

Cnsp 

0.04 

'-ma 

-0.060 

Cm£  P 

-0.080 

TABLE  I 

Linearized  Aerodynamic  Derivatives  (M=3.95) 

Circular  Airframe 

a=10°  a=20° 

-0.082  -0.111 

-0.019  -0.003 

-0.009  -0.020 

0.022  0.028 

-0.053  -0.062 

-0.016  -0.038 

-0.009  -0.022 

0.018  0.044 

0.035  0.044 

0.17  0.22 

0.04  0.05 

-0.065  -0.118 

-0.095  -0.115 

Reference  Center  of  Gravity  (C.G)  at  0.6  1 
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The  analysis  of  the  uncoupled  pitch  channel  autopilot  is  based,  on  the 
time  responses,  and  Bode  plots  of  maneuver  plane  acceleration,  body  angular 
rates  and  tail  incidence  angles,  obtained  by  OPTSYS  control  program  of  NPS. 
For  that  reason  the  differential  equations  of  the  uncoupled  pitch  channel 
autopilot  must  be  put  In  state-space  form   X=FX+6U,  which  Is  obtained  as 

follows: 

a.  Equations  of  Control  Law 

(1)  Acceleration  Filter  Equation 

1 
X,  = nz  (II.C.4-1) 

s 

+  1 

150 
or 
X,  S+  150Xt  =  150  r^  (II.C.4-2) 

and  utilizing  inverse  Laplace  Transformation  it  becomes: 
X,  =  -150Xt  ♦  150^  (II.C.4-3) 

(2)  Acceleration  Compensator  Equation 

s  + 0.143 
(-0.0387)  


0.143 


(II.C.4-4) 


(n^-X,)  s  (s  +  5) 


5 
or  by  minor  manipulation 


-0.27063  s  -0.0387 


(II.C.4-5) 


(nzc-X,)  0.2s2+s 

and  utilizing  inverse  Laplace  Transformation  it  becomes: 

36 


Y+5Y=- 1 .353 1  (n^-X,  )-0. 1 9Z5(r\Z0-X^ )  (II.C.4-6) 

Using  state-space  representation  of  a  system,  in  which  the 
forcing  function  involves  derivatives  terms  [Ref.  5:  pp  675-678]  one  yields: 

X2=  X3- 1 .353 1  (n^-X, )  (II.C.4-7) 

X3=  -5X3+6.572(nzc-X1 )  (IIX.4-8) 

Y  =  X2  (IIX.4-9) 
(3)  Rate  Compensator  Equation 

s  +  6.55 
(-  15.6) 


6.55 
Spc  = (  X2  -  X4  )  (II.C.4- 1 0) 

S  + 0.143 


0.143 


or 


-2.3816  s-  15.6 
Spc  = (  X2  -  X4  )  (II.C.4- 1 1 ) 

6.993  s  +  1 

and  using  inverse  Laplace  Transformation  it  yields: 
SpC=-0. 1 436^-340.58  x  10-3(X2-X4)-2.2308(X2-X4)  (II.C.4-12) 
b.  Actuator  Equation 

1 

Sp  = Spe  (II.C.4- 1 3) 

s+  188.4 


188.4 
and  using  inverse  Laplace  Transformation  it  yields: 
5p  =  - 1 88.4  Sp+ 1 88.4  8^  (II.C.4- 1 4) 
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c.  Equations  of  Aerodynamic  Model 

From  aerodynamic  transfer  function,  Eq.  (II.C.3-3) 


-1.3361  s- 0.159 
,   or 


6p  22.545  X  1 0"3  s2  +  3.3634  x  1 0"3  s  +  1 

22.545  x  1 0-3  q  S2  +  3.3634  x  10"3qs  +  q  =  -1.3361  SPs- 0.159  6p 

(II.C.4-15) 

and  using  inverse  Laplace  Transformation  it  yields: 

q  +  149.18610-3  q  +  44.3557  q  =-59.2637  6p  -7.0526  8p   (II.C.4- 1 6) 

Using  state-space  representation  of  a  system,  In  which  the  forcing 
function  involves  derivatives  terms  [Ref.  5:  pp  675-678]  one  obtain: 

X4  =  X5  -59.2637  6p  (II.C.4- 1 7) 

Xs  =- 1 49. 1 86  x  1 0-3  X5  -44.3557  X4  +  1 .7887  Sp  (II.C.4- 1 8) 

X4  =  q  (II.C.4-19) 

From  aerodynamic  transfer  function  Eq.(II.C.3-4) 

ITj             1 06.074  X  10"3s2-  18.82 
= ,  or 


6p  22.545  x  1 0-3  s2  +  3.3634  x  1 0"3  3  +  1 

22.545  x  1 0-3  s2  r^  +  3.3634  x  1 0*3  s  r^  +  n2=  1 06.074  x  1 0"3  s2  Sp- 1 8.828p 

(II.C.4-20) 
or  1^+1 49. 1 86  x  1 0-3  rij  +  44.3557  r^  =  4.705  Sp  -  834.775  8p 

(II.C.4-21) 

and  proceeding  as  above  it  yields: 

X*  =  X7  -  70 1 .92  x  1 0-36p  (II.C.4-22) 
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X7  =  - 1 49. 1 86  x  1 0-3  x7  -  44.3557  X6  -  1 043.366p  (II.C.4-23) 

l\  =  X6  +  4.7056p  (II.C.4-24) 

Utilizing    state-space    representation,    all    the    aforementioned 

equations  can  be  modeled  in  a  tenth  order  system  as  follows: 

X1=-150Xt  +  150X6+ 705.75  5p  (II.C.4-25) 

X2  =  X3  -  1 .353 1  njC  +  1 .353 1  X,  (II.C.4-26) 

X3  =  -5  X3  +  6.572  nzc  -  6.572  X,  (II.C.4-27) 

X4  =  X5  -  59.2637  6p  (II.C.4-28) 

X«5  =  - 1 49. 1 86  X 1 0-3  X5  -44.3557  X4  +  1 .7887  8p  (II.C.4-29) 

X*  =  X7  -  70 1 .92  x  1 0-3  Sp  (II.C.4-30) 

X7  =  - 1 49. 1 86  x  1 0-3  X7  -  44.3557  X6  -  1 043.36  Sp  (II.C.4-3 1 ) 

Sp  =  - 1 88.4  Sp  +  1 88.4  S^  (IIX.4-32) 

Spo  =  -0.143  6pc  -2.2308  X2  +  2.2308  X4  -340.58  x  10"3  X3  + 

+  460.839  x  10-3  nzc  -  460.839  x  10"3  x,  +  340.58  x  10"3  X^  - 

-20.184Sp  (II.C.4-33) 

rij  =  -(  X7  -  887. 1 24  5p  +  886.422  8^ )  (II.C.4-34) 

A  continuous,  linear,  constant  coefficient  system  of  differential 

equations,  as  above,  can  always  be  expressed  as  a  set  of  first-order  matrix 

differential  equations  X  =  FX+GU    where  U  is  the  control  input  to  the 

system.  The  output  can  be  expressed  as  a  linear  combination  of  the  state,  X, 

and  the  input  as  y=  HX+JU 

The  state  method  of  representing  a  dynamic  system  is  very  useful 
because  it  standardizes  the  information  required  into  three  matrices,  F,  G, 
and  H,  no  matter  how  complicated  the  system  is. 
In  our  case  the  matrices  are: 
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x= 


x, 

*3 

x4 

*5 
*7 

K 


x= 


x, 
x2 

Xj 

X4 

x5 
x« 

X7 

«p 

ope 


G= 


0 

-1.3531 

6.572 

0 

0 

0 

0 

0 

0.4608 

0 


U  =  n« 


F= 


-150 

0 

0 

0     0    0 

0 

0 

0 

150 

1.353 

0 

1 

0     0    0 

0 

0 

0 

0 

-6.572 

0 

-5 

0     0    0 

0 

0 

0 

0 

0 

0 

0 

0     1    0 

0 

-59.26 

0 

0 

0 

0 

0 

-44.36-0.1492  0 

0 

1.789 

0 

0 

0 

0 

0 

0     0    0 

1 

-0.7019 

0 

0 

0 

0 

0 

0     0  -44.36 

-0. 1 492 

-1043 

0 

0 

0 

0 

0 

0     0    0 

0 

-188.4 

188.4 

0 

■0.4608 

-2.231 

-0.3406 

2.231  0.3406  0 

0 

-20.18 

-0. 1 43 

0 

0 

0 

0 

0     0    0 

-1 

887.1 

-886.4 

0 

The  state-variables  of  interest  are: 
X4  =  q  :  pitch  angular  rate 

nz        :  pitch  normal  acceleration 

Sp        :  pitch  tail  incidence 

Utilizing  OPTSYS  the  tenth  order  system  was  analyzed  in  terms  of 
time  and  frequency  responses,  using  an  input  step  function  representing  the 
"1  gee  command"  and  trim  angle-of-attack  ae=0,  and  the  corresponding  plots 

were  obtained. 
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Figure  2.7  shows  the  pitch  normal  acceleration  (n3)  response  which 
has  a  0.5  seconds  time  constant,  steady-state  error  0.005  and  2.0  % 
overshoot,  and  therefore  is  according  to  the  requirements  mentioned  in 
section  II.A. 

Figures  2.8  and  2.9  show  the  required  body  angular  rate  (q)  and 
control  surface  deflection  (Sp),  to  achieve  the  acceleration  response.  These 

responses  matched  those  presented  in  [Ref.  3]. 

Figures  2.10  through  2.15  show  the  frequency  responses  of  normal 
acceleration  (r^),  angular  rate  (q),  and  tail  incidence  (5p).  The  relative 

stability  margins  are  shown  in  Table  II.  Since  the  gain  and  phase  margins,  as 
shown  in  Table  II  are  positive,  it  ensures  that  the  closed  loop  (controlled) 

system  will  be  stable. 

Figures  2.16  through  2.18  show  the  Pole-Zero  Maps  of  n2,  q  and    Sp  The 

poles  locations  are  found  to  be  at: 


St   =-159.7  + j  18.99 

(Acceleration  Filter) 

(II.C.4-35) 

Sj  =-159.7-]  18.99 

(Acceleration  Compensator) 

(II.C.4-36) 

s3  =  -8.290  *  j  8.059 

• 

(II.C.4-37) 

s4  =  -8.290  -  j  8.059 

(Angular  Rate) 

(II.C.4-38) 

%  =-3.759  + j  2.515 

(II.C.4-39) 

%  =-3.759- j  2.515 

(II.C.4-40) 

s7  =  -0. 1 4390000 

(II.C.4-41) 

Sg  =-0.000004251 

(Actuator) 

(II.C.4-42) 

%  =  -0.07459  +  j  6.66 

(Rate  Compensator) 

(II.C.4-43) 

s10  =  -0.07459  -  j  6.66 

(Normal  Acceleration) 

(II.C.4-44) 

sinse  all  poles  are  in  the  left-half  of  the  s-plane  the 

system  is  stable. 
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Figure  2.7    Pilch  Normal  Acceleration  vs  Time;  Uncoupled  Pitch  Channel; 
Continuous  Open  Loop  System. 
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Figure  2.8  Pitch  Angular  Rate  vs  Time;  Uncoupled  Pitch  Channel; 
Continuous  Open  Loop  System. 
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Figure  2.9    Pitch  Tall  Incidence  vs  Time;  Ucoupled  Pitch  Channel; 
Continuous  Open  Loop  System. 
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Figure  2.10  Pitch  Normal  Acceleration;  Uncoupled  Pitch  Channel; 
Continuous  Open  Loop  System,  Gain  vs  Frequency. 
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Figure  2.11  Pitch  Normal  Acceleration;  Uncoupled  Pitch  Channel; 
Continuous  Open  Loop  System;  Phase  vs  Frequency. 
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Figure  2.12  Pitch  Angular  Rote;  Uncoupled  Pitch  Channel; 

Continuous  Open  Loop  System;  Gain  vs  Frequency. 
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Figure  2.13  Pitch  Angular  Rote;  Uncoupled  Pitch  Channel; 

Continuous  Open  Loop  System;  Phase  vs  Frequency. 
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Figure  2.14  Pitch  Tail  Incidence;  Uncoupled  Pitch  Channel; 

Continuous  Open  Loop  System;  Gain  vs  Frequency. 
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Figure  2.15  Pilch  Toll  Incidence;  Uncoupled  Pitch  Channel; 

Continuous  Open  Loop  System;  Phase  vs  Frequency. 
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TABLE   II 

Stability  Margins  of  Uncoupled  Pitch  Channel 

q  Sp  nz 

Phase  Crossover  Frequency  (rad/sec):  56.6896  57.2 1 5 1 

Gain  Margin  (DB):                                  29.0988  7.5450 

Gain  Crossover  Frequency  (rad/sec):     5217.1016  375.9622 

Phase  Margin  (Degrees):                       87.9482  62.9508 
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Figure  2.16  Pitch  Normal  Acceleration;  Uncoupled  Pitch  Channel; 
Continuous  Open  Loop  System;  Pole-Zero  Nap. 
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Figure  2.17  Pitch  Angular  Rote;  Uncoupled  Pitch  Channel; 
Continuous  Open  Loop  System;  Pole-Zero  Map. 
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Figure  2.18  Pitch  Toll  Incidence;  Uncoupled  Pitch  Channel; 
Continuous  Open  Loop  System;  Pole-Zero  flop. 
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The  zeros  of  pitch  normal  acceleration  are  found  to  be  at: 

s,   =-150.00  (II.C.4-45) 

52  =    13.32  (II.C.4-46) 

53  =-13.32  (II.C.4-47) 
s4  =-6.550  (II.C.4-48) 

55  =-0.000008758  (II.C.4-49) 

56  =     0.000004509  (II.C.4-50) 

57  =-0.1430  (II.C.4-51) 
Sq  =  -  0.07459  +  j  6.66  (II.C.4-52) 
%  =  -  0.07459  -  j  6.66  (II.C.4-53) 

The  zeros  of  pitch  angular  rate  are  found  to  be  at: 

st   =-150.00  (II.C.4-54) 

52  =-0.119  (II.C.4-55) 

53  =  0.001734  (II.C.4-56) 
s4  =-6.550  (II.C.4-57) 
S5  =  -  0.000867  +  j  0.00 1516  (II.C.4-58) 
s6  =  -  0.000867  -  j  0.00 1516  (II.C.4-59) 
s7  =-0.1430  (II.C.4-60) 
s8  =  -  0.07459  +  j  6.66  (II.C.4-6 1 ) 
%  =  -  0.07459  -  j  6.66  (II.C.4-62) 

The  zeros  of  pitch  tail  incidence  are  found  to  be  at: 

^   =-150.00  (II.C.4-63) 

S2  =-0.000001322  (II.C.4-64) 

s3  =-6.550  (II.C.4-65) 

s4  =    0.000001189  (II.C.4-66) 
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S3  =-0.1430  (II.C.4-67) 

s6  =  -  0.07459  +  j  6.66  (II.C.4-68) 

s7  =  -  0.07459  -  j  6.66  (II.C.4-69) 

s8  =  -  0.07459  *  j  6.66  (II.C.4-70) 

S9  =  -  0.07459  -  ]  6.66  (II.C.4-7 1 ) 

5.  Design  Approach  and  Analysis  of  the  Discrete  Open  Loop  System 

The  control  of  physical  systems  with  a  digital  computer  is  becoming 
more  and  more  common.  Many  new  digital  control  applications  are  being 
stimulated  by  microprocessor  technology.  Among  the  advantages  of  digital 
logic  for  control  are  the  increased  flexibility  of  the  control  programs  and 
the  decision-making  or  logic  capability  of  digital  systems  which  can  be 
shared  with  the  control  function  to  meet  other  system  requirements. 

In  order  to  compute  the  state  X(t)  by  use  of  a  digital  computer,  the 
continuous-time  state  equation  X  =  FX+GU  must  be  converted  into  the 
discrete-time  equation  X  =  AX+BU  ,  assuming  that  the  input  vector  U 
changes  only  at  the  equally  spaced  sampling  instants  T. 

Utilizing  the  Analog-to-Digital  conversion  of  ORACLS  the  continuous 
time  matrices  F  and  G  were  converted  to  the  discrete  time  matrices  A  and 
B  using  sample  time  T=0.0125  sec  (Frequency=80  cps). 

The  computer  listing  of  the  above  conversion  is  shown  in  Appendix  B. 

The  TRANFUNC  option  of  ORACLS  was  used  to  analyze  the  dicrete  open 
loop  transfer  function  and  time  responses  and  pole-zero  maps  plots  were 
obtained  as  for  the  continuous  system. 

Figures  2.19  shows  the  discrete  time  response  of  pitch  normal 
acceleration  (r^)  which  has  a  0.5  seconds  time  constant,  steady-state  error 
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Figure  2.19  Pitch  Normal  Acceleration  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Open  Loop  System. 
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about  0.005  and  2.0S  overshoot,  and  therefore  is  according  to  the 
requirements  mentlonted  in  section  II.A. 

Figures  2.20  and  2.21  show  the  discrete  time  responses  of  pitch 
angular  rate  (q),  and  pitch  tail  incidence  (6p).  These  responses  matched 

those  presented  in  [Ref.  3].  An  input  step  function  representing  "1  geeM  was 

used  in  all  cases. 

Figures  2.22  through  2.24  show  the  discrete  Pole-Zero  Maps  of  r\z,  q, 

and  Sp.  The  poles  are  found  to  be  at: 

zt  =0.1 320  +  ]  0.03 1 94    (Acceleration  Filter)  (II.C.5- 1 ) 

z2  =  0. 1 320  -  j  0.03 1 94    (Acceleration  Compensator)  (II.C.5-2) 

23  =  0.8970  +  j  0.09067  (II.C.5-3) 

z4  =  0.8970  -  j  0.09067    (Angular  Rate)  (II.C.5-4) 

Z5  =  0.9536  +  j  0.02999  (II.C.5-5) 

Z6  =  0.9536  -  j  0.02999  (II.C.5-6) 

z7  =  0.9982  (II.C.5-7) 

z8=  1.0000  (Actuator)  (II.C.5-8) 

z9  =  0.9956  +  ]  0.08307    (Rate  Compensator)  (II.C.5-9) 

z10=  0.9956  -  ]  0.08307   (Normal  Acceleration)  (II.C.5-10) 
Since  all  the  poles  are  within  the  unit  circle  the  system  is  stable. 
The  zeros  of  the  discrete  pitch  normal  acceleration  are  found  to  be  at: 

z,   =-0.4663  (II.C.5-11) 

z2  =     0.1532  (II.C.5-12) 

Z3  =     0.8466  (II.C.5-13) 

z4  =      1.1810  (II.C.5-14) 

z5  =     0.9214  (II.C.5-15) 
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Figure  2.20  Pitch  Angular  Rate  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Open  Loop  System. 
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Figure  2.21  Pitch  Toil  Incidence  vs  Time;  Ucoupled  Pitch  Channel; 
Discrete  Open  Loop  System. 
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Figure  2.22  Pitch  Normal  Acceleration;  Uncoupled  Pitch  Channel; 
Discrete  Open  Loop  System;  Pole-Zero  Map. 
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Figure  2.23  Pitch  Angular  Rate;  Uncoupled  Pitch  Channel; 
Discrete  Open  Loop  System;  Pole-Zero  Map. 
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Figure  2.24  Pitch  Tall  Incidence;  Uncoupled  Pitch  Channel; 
Discrete  Open  Loop  System;  Pole-Zero  Map. 
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26    = 

0.9982 

(II.C.5-16) 

27    = 

1.0000 

(II.C.5-17) 

28    = 

0.9956  ♦  j  0.08307 

(II.C.5-18) 

29   = 

0.9956  -  j  0.08307 

(II.C.5-19) 

The  zeros  of  the  discrete  pitch  angular  rate  are 

found  to  be  at: 

21    =" 

■  2.2750 

(II.C.5-20) 

22  = 

0.1534 

(II.C.5-21) 

23  =" 

■0.1386 

(II.C.5-22) 

24  = 

0.9985 

(II.C.5-23) 

25  = 

0.9214 

(II.C.5-24) 

26    = 

0.9982 

(II.C.5-25) 

2?    = 

1.0000 

(II.C.5-26) 

28   = 

0.9956  +  j  0.08307 

(II.C.5-27) 

29   = 

0.9956  -  j  0.08307 

(II.C.5-28) 

The  zeros  of  the  discrete  pitch  tail  incidence  are  found  to  be  at: 

21    =" 

-  0.4658 

(II.C.5-29) 

z2  = 

0.1532 

(II.C.5-30) 

23  = 

0.9956  -  ]  0.08307 

(II.C.5-31) 

24  = 

0.9956  ♦  j  0.08307 

(II.C.5-32) 

25  = 

0.9214 

(II.C.5-33) 

26    = 

0.9982 

(II.C.5-34) 

27    = 

1.0000 

(II.C.5-35) 

28    = 

0.9956  +  J  0.08307 

(II.C.5-36) 

29    = 

0.9956  -  j  0.08307 

(II.C.5-37) 
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6.  Comparison  of  Continuous  and  Discrete  Systems 

In  this  chapter  the  classical  design  of  the  continuous  uncoupled  pitch 
channel  autopilot,  as  was  developed  in  [Ref.  3],  was  reviewed,  and  the 
corresponding  discrete  channel  was  obtained  using  state-space 
representation  of  the  system  and  the  Analog-to-Digital  Conversion  option 
of  ORACLS. 

A  comparison  of  the  time  responses  in  the  s-domain  (Figures  2.7 
through  2.9)  with  those  in  the  z-domaln  (Figures  2.19  through  2.21)  shows 
that  the  two  systems  behave  the  same  way.  In  addition  both  systems  are 
stable  as  it  can  be  seen  from  the  pole-zero  maps  of  the  continuous  system 
(Figures  2.13  through  2.15) ,  and  the  discrete  system  (Figures  2.22  through 
2.24) ,  since  their  poles  are  maping  according  to  relation  z=esT. 

The  similarity  of  the  continuous  and  discrete  systems  allow  us  from 

now  on  to  work  on  the  discrete  domain  without  any  loss  of  accurasy,  and 
using  the  computer  to  make  things  much  easier. 
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III.  MODERN  CONTROL  DESIGN  AND  ANALYSIS  OF  LONGITUDINAL 
UNCOUPLED  PITCH  CHANNEL  AUTOPILOT 


A.  GENERAL 

The  design  technique  for  the  uncoupled  pitch  channel,  as  developed  in 
[Ref.  31,  was  classical  using  a  combination  of  frequency  response  and  root- 
locus  techniques. 

In  classical  control  theory,  only  the  Input,  output,  and  error  signals  are 
considered  important;  the  analysis  and  design  of  control  systems  are 
carried  out  using  transfer  functions,  together  with  a  variety  of  graphical 
techniques  as  root-locus  plots  and  Nyquist  plots.  The  unique  characteristic 
of  classical  control  theory  is  that  it  Is  based  on  the  Input-output  relation  of 
the  system,  or  the  transfer  function. 

The  main  disadvantage  of  classical  control  theory  is  that,  generally 
speaking,  it  is  applicable  only  to  linear  time-invariant  systems  having  a 
single  input  and  a  single  output.  It  is  powerless  for  time-varying  systems, 
non-linear  systems  (except  simple  ones),  and  multiple-input-multiple- 
output  systems.  Thus  classical  techniques  (the  root-locus  and  frequency 
response  methods)  do  not  apply  to  the  design  of  optimal  and  adaptive  control 
systems,  which  are  mostly  time-varying  and/or  linear. 

The  modern  trend  in  engineering  systems  Is  toward  greater  complexity, 
due  mainly  to  the  requirements  of  complex  task  and  good  accuracy.  Complex 
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systems  may  have  multiple  inputs  and  multiple  outputs  and  may  be  time- 
varying. 

Because  of  the  necessity  of  meeting  increasingly  stringent  requirements 
on  the  performance  of  control  systems,  the  increase  in  system 
complexity,and  easy  access  to  large-scale  computers,  modern  control 
theory,  which  is  a  new  approach  to  the  analysis  and  design  of  complex 
control  systems,  has  been  developed  since  around  1960.  This  new  approach 
is  based  on  the  concept  of  state-space  representation. 

Modern  control  theory  is  contrasted  with  classical  control  theory  in  that 
the  former  is  applicable  to  multiple-input-multiple-output  systems,  which 
may  be  linear  or  non-linear,  time-invariant  or  time-varying,  while  the 
latter  is  applicable  only  to  linear  time-invariant  single-input-single-output 
systems.  Also,  modern  control  theory  is  essentially  a  time-domain 
approach,  while  classical  control  theory  is  a  complex  frequency-domain 
approach. 

System  design  in  classical  control  theory  is  based  on  trial -and-error 
procedures  which,  in  general,  will  not  yield  optimal  control  systems. 
System  design  in  modern  control  theory,  on  the  other  hand,  enables  the 
engineer  to  design  optimal  control  systems  with  respect  to  given 
performance  indexes.  In  addition,  design  in  modern  control  theory  can  be 
carried  out  for  a  class  of  inputs,  instead  of  a  specific  input  function,  such 
as  the  impulse  function,  step  function,  or  sinusoidal  function.  Also,  modern 
control  theory  enables  the  engineer  to  include  initial  conditions  in  the 
design. 
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The  equations  of  the  uncoupled  pitch  channel  autopilot  have  already  been 
In  state-space  form  X  =  AX  +  BU  ,in  the  previous  chapter  ,  and  therefore 
we  proceed  in  the  modern  control  design  of  the  discrete  state-feedback 
autopilot. 

B.  STATE  FEEDBACK  DESIGNED  AUTOPILOT 

1.  Design  Approach  and  Analysis  of  Control  Law 

One  of  the  attractive  features  of  state-space  design  methods  1s  that 
the  procedure  consists  of  two  Independent  steps.  One  step  assumes  that  we 
have  all  the  states  at  our  disposal  for  feedback  purposes.  In  general,  of 
course,  this  would  be  unrealistic  since  a  practical  engineer  would  not,  as  a 
rule,  find  it  necessary  to  purchase  this  large  number  of  sensors,  especially 
since  he  knows  that  he  would  not  need  them,  using  classical  design  methods. 
The  assumption  that  all  states  are  available  merely  allows  us  to  proceed 
with  the  first  step,  namely,  the  control  law.  The  remaining  step  is  to  design 
an  "estimator''  which  estimates  the  entire  state  vector,  given  measurements 
of  the  portion  of  the  state  provided  by  y  =  HX  +  JU. 

The  final  control  algorithm  will  consist  of  the  control  law  and  the 
estimator  combined  where  the  control-law  calculations  are  based  on  the 
estimated  states  rather  than  the  actual  states.  This  substitution  Is 
reasonable  and  the  combined  control  law  and  estimator  will  give  closed- 
loop  dynamic  characteristics  which  are  unchanged  from  those  assumed  in 
designing  the  control  law  and  the  estimator  separately.  The  dynamic  system 
we  obtain  from  the  combined  control  law  and  the  estimator  is  called  the 
control ler.The  first  step  is  to  get  a  good  control  law. 
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The  control  law  is  simply  the  feedback  of  a  linear  combination  of  all 
the  states,  that  is, 

U  =  -KX  (III.B.1-0 

The  characteristic  equation  of  the  controlled  (closed  loop)  system  is 
det[  zl-A+BK  ]  =  O  (III.B.1-2) 

The  control  law  design  then  consists  of  picking  the  elements  of  K  so 
that  the  roots  of  equation  III.B.1-2  are  in  desirable  locations. 

A  program  logic  for  application  of  Ackermann's  formula  to  compute 
the  control  law  is  given  in  Appendix  C.  The  program  reads  in  the  sample  time 
T,  the  discrete  matrices  A  ,  B  and  the  desired  pole  locations  in  the  s-plane 
and  gives  as  output  the  control  law  matrix  K. 

The  program  used  in  this  work  was  written  in  WATFIV  computer 
language  and  is  given  in  Appendix  D. 

In  section  III.C.5  utilizing  Analog-to-Dlgltal  conversion  the  discrete 
uncoupled  pitch  channel  matrices  A  and  B  were  obtained.  Also  the  analysis 
of  the  continuous  open  loop  system,  obtained  from  OPTSYS,  showed  the 
eigenvalues  of  the  classical  design  to  be  at: 

s,  =  - 1 59.724  +  j  1 8.992 1  (III.B.  1  -3) 

Sj  =  - 1 59.724  -  j  1 8.992 1  (III.B.  1  -4) 

S3  =  -8.29048  +  j  8.05932  (III.B.  1  -5) 

s4  =  -8.29048  -  j  8.05932  (III.B.  1  -6) 

S5  =  -3.75925  +  j  2.51 463  (III.B.  1  -7) 

%  =  -3.75925  -  j  2:5 1 463  (III.B.  1  -8) 

s7  =-0.14393000  (III.B.  1-9) 

Se  =-0.00000425  (III.B.  1-10) 
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S9  =  -0.074593  ♦  J  6.65959  (III.B.  1-11) 

s10=  -0.074593  -  j  6.65959  (III.B.  1-12) 

Utilizing     Ackermann's  formula  program  (Appendix  D)  and  using 
sample  time  T=  0.0125  sec  (80  cps),  the  discrete  matrices  A,B  (Appendix  B) 
and  desired  poles  locations  in  s-plane  the  continuous  open  loop  eigenvalues 
(equations  III.B.  1-3  through  III.B.  1-12)  the  control  law  was  found  to  be: 
K  =  1-0.0188  -0.0279         -0.0042         0.0290  0.0041 

0.0  0.0001  1.7732  -2.0162         0.0131] 

(III.B.  1-13) 

Utilizing  Transfer  Function  Analysis  from  TRANFUNC  option  of 
ORACLS,  and  using  the  discrete  matrices  A,B  and  K  the  state-feedback 
autopilot  was  designed  and  analysed  according  to  step  responses  and  Pole- 
Zero  Maps. 

2.  Performance  of  System 

A  comparison  of  the  discrete  open  loop  eigenvalues  (detl  zl-A  ]=0), 
obtained  using  the  Compensator  Transfer  Function  Analysis  from  TRANFUNC 
option  of  ORACLS: 


zt  =0.131993  + j  0.0319377 

Z2  =0.131993- j  0.0319377 

Z3  =  0.896987  +  ]  0.0906708 

z4  =  0.896987  -  J  0.0906708 

Z5  =  0.953625  +  j  0.0299850 

Z6  =  0.953625  -  j  0.0299850 

Z7  =0.998202 

z8  =  1.000000 


(III.B.2-1) 
(III.B.2-2) 
(III.B.2-3) 
(III.B.2-4) 
(III.B.2-5) 
(III.B.2-6) 
(III.B.2-7) 
(III.B.2-8) 
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Zg  =  0.995608  ♦  j  0.0830713  (III.B.2-9) 

Z10  =  0.995608  -  j  0.08307 1 3  (III.B.2- 1 0) 

with  the  closed   loop   regulator   eigenvalues  (det[   zl-A+BX    1=0) 

obtained  using  the  Transfer  Function  Analysis  from  TRANFUNC  option  of 

ORACLS: 

Zj  =0.134462  +  ]  0.0320974 
Z2  =0.134462-]  0.0320974 
Z3  =  0.896990  +  j  0.0906457 
Z4  =  0.896990  -  j  0.0906457 
Z5  =  0.953624  +  j  0.0299994 
Z6  =  0.953624  -  j  0.0299994 
Z7  =  0.998202 
Z8  =  1.000000 

Z9  =0.995608  +  ]  0.0830713 
Z10  =  0.995608 -j  0.0830713 

shows  no  significant  differences. 

Figures  3.1  through  3.3  show  the  discrete  closed  loop  time  responses 
of  pitch  normal  acceleration  (nz),  pitch  angular  rate  (q),  and  pitch  tail 
incidence  (Sp).  An  input  step  function  representing  "1  gee"  was  used.  A 
comparison  of  the  responses  mentioned  above,  with  the  corresponding  time 
responses  of  the  discrete  open  loop  system  figures  2.19  through  2.21  shows 
that  they  are  identical. 

Figures  3.4  through  3.6  show  the  Pole-Zero  Map  plots  of  nz,  q,  and  5p. 

All  the  poles,given  by  equations  III.B.2- 1 1  through  III.B.2-20,  are  inside  the 
unit  circle  and  the  system  is  stable.  There  are  no  real  zeros  in  the  system. 


(III.B.2-11) 
(III.B.2- 12) 
(III.B.2- 13) 
(III.B.2- 14) 
(III.B.2- 15) 
(III.B.2- 16) 
(III.B.2- 17) 
(III.B.2- 18) 
(III.B.2- 19) 
(III.B.2-20) 
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Figure  3.1     Pitch  Normal  Acceleration  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  State-Feedback  Autopilot. 
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Figure  3.2 


Pitch  Angular  Rate  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  State-Feedback  Autopilot. 
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Figure  3.3 


Pitch  Tall  Incidence  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  State-Feedback  Autopilot. 
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Figure  3.4    Pitch  Normal  Acceleration  ;  Uncoupled  Pitch  Channel; 
Discrete  State-Feedback  Autopilot;  Pole-Zero  Map. 
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Figure  3.5    Pitch  Angular  Rote;  Uncoupled  Pitch  Channel; 

Discrete  State-Feedback  Autopilot;  Pole-Zero  Map. 
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Figure  3.6    Pitch  Toil  Incidence;  Uncoupled  Pitch  Channel; 

Discrete  State-Feedback  Autopilot;  Pole-Zero  Map. 
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3.  Simplified  State-Feedback  Designed  Autopilot 

In  the  previous  section  the  state-feedback  autopilot  was  designed 
using  the  feedback  gain  vector  (equation  III.B.1-13),  and  therefore  there 
were  ten  (10)  returning  gain  loops  making  the  system  more  complex.  In 
order  to  simplify  the  system,  zero  gains  were  Introdused,  where  those  gains 
were  very  small  in  the  state  feedback  gain  vector,and  In  doing  so  returning 
gain  loops  were  eliminated  and  the  system  became  simpler. 
The  resulted  state-feedback  gain  vector  Is: 
K  =  [-0.0188  -0.0279         0.0000  0.0290        0.0000 

0.0000  0.0000         1.7732  -2.0162         0.0131] 

(IILB.3-1) 
Following  the  same  procedure,  as  In  the  case  of  State-Feedback 
Designed  Autopilot,  the  step  responses  and  Pole-Zero  Map  plots  were 
obtained  and  analysed. 

The  pole  locations  of  the  simplified  state-feedback  autopilot  are 
found  to  be  at: 

2t  =  0. 1 34458  ♦  j  0.032 1 064  (III.B.3-2) 

z2  =  0. 1 34458  -  ]  0.032 1 064  (III.B.3-3) 

Z3  =  0.896805  +  j  0.09 1 1 1 70  (III.B.3-4) 

z4  =  0.896805  -  j  0.09 1 1 1 70  (III.B.3-5) 

Z5  =  0.953644  +  ]  0.0290202  (III.B.3-6) 

z6  =  0.953644  -  j  0.0290202  (III.B.3-7) 

Z7  =  0.998202  (III.B.3-8) 

Z8=  1.000000  (III.B.3-9) 

Z9  =  0.995608  +  j  0.08307 1 3  (III.B.3- 1 0) 
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Z10=  0.995608  -  ]  0.08307 1 3  (III.B.3- 1 1 ) 

A  comparison  of  the  above  pole  locations  with  those  of  the  state- 
feedback  designed  autopilot  (equations  III.B.2-1 1  through  III.B.2-20)  show  a 
very  small  shifting  of  the  poles  of  the  simplified  state-feedback  autopilot 
to  the  left. 

0 

Figure  3.7  shows  the  time  response  of  the  pitch  normal  acceleration 
(r^)  for  the  simplified  state-feedback  autopilot.  An  Input  step  function 

representing  "1  gee"  was  used.  The  response  shows  that  the  steady-state 
error  is  greater  than  in  the  case  of  state-feedback  autopilot  (Fig.  3.1),  but 
is  whithin  the  requirements. 

Figures  3.8  and  3.9  show  the  time  responses  of  pitch  angular  rate  (q) 
and  pitch  tail  incidence  (Sp)  for  the  simplified  state-feedback  autopilot.  A 

comparison  with  the  corresponding  time  responses  for  the  state-feedback 
autopilot  (Fig.  3.2  and  Fig.  3.3)  shows  no  significant  differences. 

Figures  3.10  through  3.12  show  the  Pole-Zero  Map  plots  of  n2,  q,  8p. 

All  the  poles,  given  by  equations  III.B.3-2  through  III.B.3- 1 1,  are  inside  the 
unit  circle  and  the  system  is  stable.  There  are  no  real  zeros  in  the  system. 

C.  ESTIMATOR  DESIGNED  AUTOPILOT 

1.  Design  Approach  and  Analysis  of  Discrete  Closed  Loop  Estimator 
The  control  law  designed  in  the  last  section  III.B  assumed  that  all 
states  were  available  for  feedback.  Since,  typically,  not  all  states  are 
measured,  an  estimator  must  be  design  in  order  to  estimate  all  the  states, 
given  measurements  of  a  portion  of  them.  If  the  state  is  X,  then  the 
estimate  is  X  and  the  idea  is  to  let  U  =  -LX,  replacing  the  true  states  by 
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Figure  3.7    Pitch  Normal  Acceleration  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Simplified  State-Feedback  Autopilot. 
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Figure  3.8    Pitch  Angular  Rote  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Simplified  Stote-Feedbock  Autopilot. 
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Figure  3.9    Pitch  Tall  Incidence  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Simplified  State-Feedback  Autopilot. 
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Figure  3.10  Pilch  Normal  Acceleration  ;  Uncoupled  Pitch  Channel; 

Discrete  Simplified  State-Feedback  Autopilot;  Pole-Zero  hap. 
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Figure  3.11  Pitch  Angular  Rate;  Uncoupled  Pitch  Channel; 

Discrete  Simplified  State-Feedback  Autopilot;  Pole-Zero  Map. 


84 


e 

— • 

o 

o- 

o 

CM 

o 

x    1          ! 

x   X 

x 

© 

N 

X                       | 

xx 
x  x 

1 

i               i 

1 

CO 

i               i 

1 

S3 

i               ; 

1 

o 

1 

■    ' 

.  _ 

■ 1 r             -    1 

-1.00  -0.75 


-0.50 


-0.25 


0.00 


0.25 


0.50 


0.75 


1.00 


Figure  3.12  Pitch  Tail  Incidence;  Uncoupled  Pitch  Channel; 

Discrete  Simplified  State-Feedback  Autopilot;  Pole-Zero  Map. 
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their  estimates  In  the  control  law.  Figure  3.13  shows  an  open  loop  estimator 
and  Figure  3.14  shows  a  closed  loop  estimator. 

The  characteristic  equation  of  the  closed  loop  estimator  is: 

detl  zl-A+LH  ]  =  0  (III.C.1-1) 

To  select  L,  we  take  the  same  approach  as  we  did  when  designing  the 
control  law  K. 

If  we  take  the  transpose  of  A-LH,  we  get  AT-HTLT,  which  Is  the  same 
form  as  the  system  matrix  A-BK  of  the  control  problem.  Therefore  If  we 
substitute  AT  for  A,  HT  for  B,  and  LT  for  K,  we  can  use  the  control  design 
results. 

Then  utilizing  Ackermann's  formula  program  (Appendix  D),  and  using 
sample  time  T=  0.0125  sec,  the  matrices  AT  for  A,  HT  for  B,  where 

HT  =  (0         0         0  1  0         0         0         0         0         01 

(III.C.  1-2) 

assuming  that  we  can  only  measure  the  state  X4  =  q,  and  taking  the 

desired  pole  locations  of  the  estimator  to  be  slightly  to  the  left  from  the 
original  open  loop  eigenvalues  (equations  III.B.1-3  through  III.B.1-12),  as 

shown  below: 

s,  =  - 1 59.729  +  j  1 8.992 1  (III.C.  1  -3) 

52  =  - 1 59.729  -  j  1 8.992 1  (III.C.  1  -4) 

53  =  -8.29548  +  j  8.05932  (III.C.  1  -5) 
s4  =  -8.29548  -  ]  8.05932  (III.C.  1-6) 
S5  =  -3.76425  +  j  2.5 1 463  (III.C.  1  -7) 
s6  =  -3.76425  -  j  2.5 1 463  (III.C.  1  -8) 
S7  =-0.148930  (III.C.  1-9) 
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Figure  3.13  Open  Loop  Estimator. 
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Figure  3.14  Closed  Loop  Estimator. 
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s8  =  -0.005004 

Sc»  =  -0.079593  ♦  ]  6.65959 

s10  =  -0.079593  -  j  6.65959 
we  obtain  the  estimator  gain  vector : 
U=[    -1.8784         0.0200  -0.0972 

-0.0015         0.3860  0.0720 


(III.C.1-10) 
(III.C.1-11) 
(III.C.1-12) 


0.0219  0.0041 

0.0007  -0.3370  ] 

(III.C.1-13) 

Utilizing  Compensator  Transfer  Function  Analysis  from  TRANFUNC 
option  of  ORACLS  and  using  the  matrices  A,  H  and  L  the  closed  loop 
estimator  autopilot  was  designed  and  analyzed  according  to  step  responses 
and  Pole-Zero  Maps. 

2.  Performance  of  System 

The  closed  loop  estimator  eigenvalues,  detl  zl-A+LH  ]  =  0,  obtained 
from  Compensator  Transfer  Function  Analysis  are: 
z]   =0.140601 


z2  =0.101457 

z3  =0.897018  + 

z4  =  0.897018- 

z5  =0.953599  + 

Z6  =0.953599- 

Z7  =0.995619  + 

Z8  =  0.995619- 

Zg  =0.998203 
z10  =  0.999998 


0.0908305 
0.0908305 
0.0299821 
0.0299821 
0.0830597 
0.0830597 


(III.C.2-1) 
(III.C.2-2) 
(III.C.2-3) 
(III.C.2-4) 
(III.C.2-5) 
(III.C.2-6) 
(III.C.2-7) 
(III.C.2-8) 
(III.C.2-9) 
(III.C.2-10) 
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A  comparison  of  the  above  eigenvalues  with  those  of  the  closed  loop 
regulator  eigenvalues  det[  zl-A+BK  ]=0,  (equations  III.B.2-11  through 
III.B.2-20),  shows  an  acceptable  difference  since  we  estimated  the  entire 
state  vector  from  the  measurement  of  Just  one  state,  namely  the  pitch 
angular  rate  (q),  and  we  moved  the  desired  pole  locations  of  the  estimator 
to  the  left  of  the  original  open  loop  eigenvalues,  in  order  to  get  a  faster 
response. 

Figures  3.15  through  3.17  show  the  discrete  closed  loop  estimator 
time  responses  of  pitch  normal  acceleration  (nz),  pitch  angular  rate  (q),  and 
pitch  tail  incidence  (6p).  An  input  step  function  representing  "1  gee"  was 

used.  A  comparison  of  the  responses  mentioned  above  with  the 
corresponding  time  responses  of  the  state-feedback  autopilot  Figures  3.1 
through  3.3  shows  no  differences.  This  Is  due  to  the  fact  that  the  same 
(zero)  initial  conditions  was  used  for  both  cases. 

Figures  3.18  through  3.20  show  the  Pole-Zero  Map  plots  of  nz,  q  and  Sp 

for  both  the  state-feedback  and  the  estimator  designed  autopilots.  All  the 

poles,  given  by  equations  III.C.2-1  through  III.C.2-10,are  Inside  the  unit 

circle  and  the  system  is  stable. 

The  zeros  of  the  estimator's  pitch  normal  acceleration  are  found  to  be 

at: 

z1  =  0.1053  (III.C.2-11) 

z2  =  0.8970  +  j  0.09080  (III.C.2-1 2) 

Z3  =  0.8970  -  j  0.09080  (III.C.2-1 3) 

Z4  =  0.1413  (III.C.2-1 4) 

Z5  =  0.9956  ♦  j  0.08307  (III.C.2-1 5) 
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Figure  3.15  Pitch  Normal  Acceleration  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Estimator  Autopilot. 
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Figure  3.16  Pitch  Angular  Rate  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Estimator  Autopilot. 
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Figure  3.17  Pitch  Tail  Incidence  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Estimator  Autopilot. 
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Figure  3.18  Pitch  Normal  Acceleration  ;  Uncoupled  Pitch  Channel;  Discrete 
State  Feedback  and  Estimator  Autopilots;  Pole-Zero  Map. 
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Figure  3.19  Pitch  Angular  Rote;  Uncoupled  Pitch  Channel;  Discrete  State 
Feedback  and  Estimator  Autopilots;  Pole-Zero  Map. 
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Figure  3.20  Pitch  Toil  Incidence;  Uncoupled  Pitch  Channel;  Discrete 
State  Feedback  and  Estimator  Autopilots;  Pole-Zero  Map. 
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26    = 

0.9956  -  ]  0.08307 

2?    = 

0.8468 

28    = 

1.0000 

29    = 

0.9982 

The  zeros  of  the  estimator's  pitch 

21  = 

0.9956  +  ]  0.08307 

22  = 

0.9956  -  )  0.08307 

23  = 

0.8970  + J  0.09081 

24  = 

0.8970 -]  0.09081 

25  = 

1.0000 

26   = 

0.1056 

2?   = 

0.9982 

28   = 

0.9536  ♦  j  0.02998 

29   = 

0.9536  -  ]  0.02998 

The  zeros  of  the  estimator's  pitch 

21  = 

0.9958  +  j  0.08305 

22  = 

0.9958  -  j  0.08305 

23  = 

0.9536  ♦  j  0.02997 

24  = 

0.9536  -  j  0.02997 

25  = 

0.1413 

26   = 

1.0000 

27  = 

0.9982 

28  = 

0.8970  + j  0.09081 

29  = 

0.8970- j  0.09081 

angular  rate  are  found 


III.C.2-16) 
III.C.2-17) 
III.C.2-18) 
III.C.2-19) 

0  be  at: 

III.C.2-20) 

III.C.2-21) 

III.C.2-22) 

III.C.2-23) 

III.C.2-24) 

III.C.2-25) 

III.C.2-26) 

III.C.2-27) 

III.C.2-28) 

tail  incidence  are  found  to  be  at: 

:iII.C.2-29) 

III.C.2-30) 

:iII.C.2-31) 
III.C.2-32) 
III.C.2-33) 
:iII.C.2-34) 
;iII.C.2-35) 
:iII.C.2-36) 
[IILC.2-37) 
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3.  Simplified  Estimator  Designed  Autopilot 

For  the  same  reason,  as  In  the  case  of  simplified  state-feedback 
autopilot,  in  order  to  eliminate  returning  loops  in  the  design  of  simplified 
estimator  autopilot,  zero  gains  were  Introduced  where  the  gains  of  the 
estimator  (equation  III.C.1-13)  were  very  small.  Then  the  resulted  gain 
vector  is: 

Us  I    -1.8784         0.0200  -0.0972         0.0219  0.0000 

0.0000        0.3860  0.0720  0.0000  -0.3370] 

(III.C.3-1) 
Following  the  same  procedure,  as  in  the  case  of  the  Estimator 
Designed  Autopilot,  the  time  responses  and  Pole-Zero  Map  plots  were 
obtained  and  analyzed. 

The  pole  locations  of  the  simplified  estimator  autopilot  are  found  to 
be  at: 


z,  =0.140707 

z2  =0.101090 

z3  =0.896885  + 

Z4  =  0.896885- 

z5  =0.953882  + 

z6  =0.953882- 

Z7  =  0.995605  + 

Zq  =  0.995605- 

z9  =0.998202 
Z10  =  0.999985 


0.0938778 
0.0938778 
0.0288601 
0.0288601 
0.0830655 
0.0830655 


(IILC.3-2) 

(III.C.3-3) 

(III.C.3-4) 

(III.C.3-5) 

(III.C.3-6) 

(III.C.3-7) 

(III.C.3-8) 

(III.C.3-9) 

(III.C.3-10) 

(III.C.3-1 1) 


97 


A  comparison  of  the  above  pole  locations  with  those  of  the  estimator 
designed  autopilot  (equations  III.C.2-1  through  III.C.2-10)  shows  no 
significant  differences. 

Figures  3.21  through  3.23  show  the  simplified  closed  loop  estimator 
time  responses  of  pitch  normal  acceleration  (nz),  pitch  angular  rate  (q),  and 
pitch  tail  incidence  (Sp).  An  Input  step  function  representing  M1  gee"  was 
used.  A  comparison  of  the  above  responses  with  the  corresponding  time 
responses  of  the  estimator  designed  autopilot  (Figures  3.15  through  3.17) 

shows  no  significant  differences. 

Figures  3.24  through  3.26  show  the  Pole-Zero  Map  plots  of  nr,  q  and  5p 

for  the  simplified  estimator  designed  autopilot.  All  the  poles,  given  by 
equations  III.C.3-2  through  III.C.3- 1 1  #are  inside  the  unit  circle  and  the 
system  is  stable.  There  are  no  significant  differences  In  the  zero's 
locations  from  those  of  the  estimator  autopilot  (equations  III.C.2-1 1 
through  III.C.2-37). 
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Figure  3.21  Pitch  Normal  Acceleration  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Simplified  Estimator  Autopilots. 
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Figure  3.22  Pitch  Angular  Rate  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Simplified  Estimator  Autopilot. 
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Figure  3.23  Pitch  Tail  Incidence  vs  Time;  Uncoupled  Pitch  Channel; 
Discrete  Simplified  Estimator  Autopilot. 
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Figure  3.24  Pitch  Normal  Acceleration  ;  Uncoupled  Pitch  Channel; 

Discrete  Simplified  Estimator  Autopilot;  Pole-Zero  Map. 
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Figure  3.25  Pitch  Angular  Rote;  Uncoupled  Pitch  Channel; 

Discrete  Simplified  Estimator  Autopilot;  Pole-Zero  Map. 
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Figure  3.26  Pitch  Toil  Incidence;  Uncoupled  Pitch  Channel; 

Discrete  Simplified  Estimator  Autopilot;  Pole-Zero  Nap. 
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IV.  DISCRETE  COUPLED  PITCH  AND  ROLL  CHANNEL  AUTOPILOTS 

A.  GENERAL 

The  control  design  procedures  described  in  chapter  II  and  III  were  applied 

to  systems  with  one  input  and  one  output.  The  transfer-function  approach  in 
[Ref.  3]  is  fundamentally  limited  to  single  input/output  systems  while  the 
state-space  methods  of  chapter  III  were  limited  to  single  input/output  to 
simplify  the  procedures.  In  fact,  If  one  tries  to  apply  the  "pole-placement" 
approach  of  chapter  III  to  a  multivariate  (more  than  one  input  or  output) 
system,  all  the  procedures  work  but  the  gains,  K  or  L,  are  not  uniquely 
determined  by  the  resulting  equations.  Until  a  design  approach  is  available 
which  intelligently  uses  this  extra  freedom,  the  pole-placement  concept  for 
estimator  designs  for  systems  with  more  than  one  output  or  for  controller 
designs  for  systems  with  more  than  one  input  has  limited  value. 

In  this  chapter  the  pole  placement  and  robustness  design  program 
(POPLAR),  [Ref.  8]  was  used  as  a  tool  in  the  design  of  multivariable  systems. 
This  program  is  designed  to  employ  singular  value  analysis  and  the  use  of  an 
optimization  to  aid  in  pole  placement  control  system  design  for  linear 
multivariable  systems.  Robustness  of  the  system  is  also  considered  by 
establishing  singular  value  levels  which  correspond  to  multiloop  gain  and 
phase  margins  determined  from  the  universal  gain  phase  diagram  developed 
by  Newsom  and  Mukhapadhyay  at  NASA  Langley. 
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B.  DESIGN  APPROACH  AND  ANALYSIS  OF  THE  DISCRETE  STATE 
FEEDBACK  AUTOPILOT 

The  first  step  In  any  multivariate  design  should  be  an  attempt  to 

find  an  approximate  model  consisting  of  two  or  more  single  input/output 

model  or  else  decouple  the  control  law  matrix  K  and  the  estimator  law 

matrix    L.  This  step  will  give  better  physical  insight  into  the  important 

feedback    variables    and    may    lead    to    a    plant    description    which    is 

substantially  simpler  for  design  purposes  and  yet  yields  no  significant 

degradation  from  an  analysis  based  on  the  full  multivariate  system. 

Combining  the  pitch  and  roll  channel  autopilot,  the  roll  channel  taken 
from  [Ref.  7],  a  seventeenth  (17th)  order  system  is  obtained  and  is  shown  in 
Appendix  E. 

Utilizing  Transfer  Function  Analysis  from  TRANFUNC  option  of  ORACLS 
and  using  the  combined  pitch  and  roll  matrices  A,  B  and  K  the  discrete 
closed  loop  system  was  obtained  and  analyzed  according  to  time  responses 
and  Pole-Zero  Map  plots. 

C.  PERFORMANCE  OF  SYSTEM 

The     discrete  state-feedback  eigenvalues,  det[  zl-A+BKl=0,  of   the 

coupled  system,  obtained  from  Transfer  Function  Analysis     option  of 

ORACLS,  are: 

z1  =  0.134462  + j  0.0320974  (IV.C-1) 

Z2  =  0.134462  -  j  0.0320974  (IV.C-2) 

z3  =  0.896990  ♦  j  0.0906457  (IV.C-3) 

Z4  =  0.896990  -  j  0.0906457  (IV.C-4) 

Z5  =  0.953624  ♦  j  0.0299994  (IV.C-5) 
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26  =  0.953624  -  j  0.0299994  (IV.C-6) 

27  =  0.998202  (IV.C-7) 

28  =  1.000000  (IV.C-8) 

29  =  0.995608  +  j  0.0830713  (IV.C-9) 
210=  0.995608  -  j  0.08307 1 3  (IV.C- 1 0) 
2,,=  0.0993086  (IV.C- 11) 
212=  0.836553  +  ]  0.306 1 07  (IV.C- 1 2) 
2,3=  0.836553  -  j  0.306 1 07  (IV.C- 1 3) 
214=  0.967788  +  j  0.0334452  (IV.C- 14) 
2,5=  0.967788  -  j  0.0334452  (IV.C- 15) 
2,6=  0.895730  (IV.C- 16) 
2,7=  0.938133  (IV.C- 17) 

Since  all  the  poles  are  inside  the  unit  circle  the  coupled  system  is 
stable. 

Figures  4.1  through  4.3  show  the  closed  loop  time  responses  of  pitch 
normal  acceleration  (n2),  pitch  angular  rate  (q)  and  pitch  tail  incidence  (6p). 

An  input  step  function  representing  "1  gee"  was  used.  All  the  responses  are 
identical  to  those  obtained  for  the  uncoupled  pitch  channel  and  therethore 
they  meet  the  desired  requirements. 

Figures   4.4  through   4.6   show   the   Pole-Zero   Map   of   r\z,    q   and   6p 

respectively.  The  poles  of  the  system  are  given  by  equations  IV.C- 1  through 
IV.C- 17  and  they  are  the  same  as  the  uncoupled  closed  loop  poles  (equations 
III.B.2-1 1  through  III.B.2-20).  There  are  no  real  2eros  in  the  state  feedback 
coupled  system. 


107 


LEGEND 
a    PITCH  NORMAL  ACCELERATION. NZ 


1.5 


0.0 


0.5 


1.0 


2.0 


2.5 


3.0 


TIME  -  SEC 


Figure  4.1 


Pitch  Normal  Acceleration  vs  Time;  Coupled  Pitch  and  Roll 
Channels;  Discrete  State-Feedback  Autopilot. 
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Figure  4.2 


Pitch  Angular  Rate  vs  Time,  Coupled  Pitch  and  Roll  Channels; 
Discrete  State-Feedback  Autopilot. 
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Figure  4.3    Pitch  Tail  Incidence  vs  Time;  Coupled  Pitch  and  Roll  Channels; 
Discrete  State-Feedback  Autopilot. 
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Figure  4.4    Pitch  Normal  Acceleration  ;  Coupled  Pitch  and  Roll  Channels; 
Discrete  State-Feedback  Autopilot;  Pole-Zero  Map. 
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Figure  4.5    Pitch  Angular  Rote;  Coupled  Pitch  ond  Roll  Channels; 
Discrete  State-Feedback  Autopilot;  Pole-Zero  Map. 
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Figure  4.6    Pitch  Tail  Incidence;  Coupled  Pitch  and  Roll  Channels; 
Discrete  State-Feedback  Autopilot;  Pole-Zero  Map. 
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D.  ROBUSTNESS  OF  SYSTEM 

With  the  rising  Interest  In  multivariate  control  theory  brought  on  by 
increasingly  complex  systems  the  need  has  arisen  to  develop  design  methods 
that  will  allow  the  designer  to  specify  system  performance  while  at  the 
same  time  ensuring  relatively  high  stability  margins  or  robustness.  In  the 
single-input-single-output  (SISO)  case  the  designer  has  had  the  tools  to  do 
these  tradeoffs  In  the  form  of  Nyqulst,  Bode  and  root-locus  plots.  In  the 
multi-input-multi-output  (MIMO)  case  the  classical  methods  are  not  totally 
appropriate. 

With  the  Increased  Interest  in  MIMO  systems  numerous  methods  of  design 
have  been  employed  to  obtain  suitable  system  performance  and  robustness 
with  varying  degrees  of  success.  One  primary  method  of  design  is  to  keep 
the  plant  as  decoupled  as  possible  throughout  the  design  so  that  each 
Individual  element  may  be  controlled  Independently  and  designed  essentially 
as  a  single  loop  system.  In  another  procedure  the  multiloop  system  is 
modified  Into  a  system  that  has  diagonal  elements  that  are  much  larger  than 
any  off-diagonal  elements.  This  diagonally  dominant  system  Is  then  in  a 
form  where  conventional  Nyquist  type  techniques  can  be  employed  in  the 
analysis.  A  third  common  MIMO  design  method  Is  that  of  the  Linear  Quadratic 
(LQ)  method.  This  method  uses  a  quadratic  cost  functional  and 
optimamization  principles  to  allow  the  designer  to  design  for  various 
performance  levels  by  adjusting  the  matrix  weighting  terms  used  In  the 
cost  function.  The  major  difficulty  with  all  of  the  above  methods  is  that 
they  are  not  necessarily  robust.  This  is  especially  true  for  cross-coupling 
terms  between  loops. 
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A  simple  Interpretation  of  robustness  Is  the  ability  of  a  system  to 
tolerate  design  pertubatlons.  These  pertubations  could  be  in  the  form  of 
actuator  failures,  plant  parameter  uncertainty,  unmodeled  dynamic  or 
nonlinear  terms,  or  any  one  of  many  other  pertubations  to  the  nominal 
design  of  the  system. 

The  analysis  of  the  robustness  of  the  system  was  obtained  using  the  Pole 
Placement  and  Robustness  Design  Program  (POPLAR)  developed  In  [Ref.  81. 

Utilizing  the  above  program  and  using  the  data  file  shown  In  Appendix  F 
for  the  state  feedback  autopilot,  and  the  data  file  shown  in  Appendix  G  for 
the  estimator  autopilot,  the  system  found  to  be  robust. 

Figures  4.7  and  4.8  show  the  Minimum  Additive  Output  Singular  Value 
(SVADMO)  vs  frequency  for  the  state  feedback  and  estimator  designed 
autopilots  respectively.  It  Is  noted  that  for  very  low  frequencies  the  values 
of  SVADMO  are  above  0.80680  and  0.83976  for  each  of  the  designs. 

Figures  4.9  and  4.10  show  the  Minimum  Additive  Input  Singular  Value 
(MIN.ADD.IN.SV)  vs  frequency  for  the  state  feedback  and  estimator  designed 
autopilots  respectively. 

In  the  case  of  state  feedback  autopilot  the  ordered  complex  eigenvalues 
(Input),  and  the  ordered  computed  eigenvalues  (output)  are: 

Ordered  input  Ordered  output 

-0.01000  + j  1.00000  0.09931 

0.09928  0.13446- j  0.03210 

0.13199  -  j  0.03194  0.13446  +  J  0.03210 

0. 1 3 1 99  +  j  0.03 1 94  0.83653  -  J  0.306 1 0 

0.83522  -  J  0.30974  0.83653  +  ]  0.306 1 0 
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Figure  4.7 


SVADMO  vs  Frequency;  Coupled  Pitch  and  Roll  Channels; 
Discrete  State-Feedback  Autopilot. 
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Figure  48 


SVADMO  vs  Frequency;  Coupled  Pitch  and  Roll  Channels; 
Discrete  Estimator  Autopilot. 
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Figure  4.9 


MIN.ADD.IN.SV  vs  Frequency;  Coupled  Pitch  and  Roll  Channels; 
Discrete  State-Feedback  Autopilot. 
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Discrete  Estimator  Autopilot. 
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0.83522  +  , 

|  0.30974 

0.89380 

0.89699  -  . 

|  0.09067 

0.89699  +  , 

|  0.09067 

0.95363  -  . 

|  0.02998 

0.95363  + . 

|  0.02998 

0.96909  -  , 

|  0.03286 

0.96909  +  . 

|  0.03286 

0.99561  -. 

|  0.08307 

0.99561  *t 

|  0.08307 

0.99820 

1.00000 

0.89571 

0.89698  -  , 

|  0.09064 

0.89698  +  , 

|  0.09064 

0.93813 

0.95360  -  , 

|  0.03001 

0.95360  +  , 

j  0.03001 

0.96776  - 

1  0.03345 

0.96776  +  . 

(  0.03345 

0.99559  - 

)  0.08307 

0.99559  + 

j  0.08307 

0.99820 

1.00000 

In  the  case  of  the  estimator  autopilot  the  ordered  complex  eigenvalues 
(Input)  are  the  same,  as  In  the  state  feedback  autopilot,  and  the  ordered 
computed  eigenvalues  (output)  are: 

Ordered  output  eigenvalues 

0.07717 

0.09855 

0.25923 

0.69508 

0.83707  +  ]  0.30140 

0.83707 -]  0.30140 

0.89849 

0.90064 

0.93996 
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0.96207 

0.96562  +  ]  0.03401 
0.96562- j  0.03401 
0.99559  +  j  0.08307 
0.99559  -  j  0.08307 
0.99820 
1.00000 
1.06596 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  goal  of  this  thesis  was  to  design  and  analyze  the  discrete 
longitudinal  autopilots  for  BTT  and  STT  missiles. 

The  following  are  the  principal  conclusions  based  on  this  work. 

1.  The  continuous  and  discrete  longitudinal  autopilots  were  proved  to 
have  identical  performances. 

2.  The  State-Feedback  and  Estimator  autopilots  were  introduced  as 
additional  dynamic  designs  in  order  to  implement  control  into  the  original 
system.  Both  performances  were  found  to  meet  the  desired  requirements. 

3.  The  Simplified  State-Feedback  and  Estimator  autopilots  reduced  the 
complexity  of  the  system,  by  eliminating  returning  gain  loops,  and  without 
any  significant  effects  on  the  system's  performance. 

4.  The  performance  of  the  Coupled  Pitch  and  Roll  channel  autopilots  was 
found  to  be  satisfactory  and  also  proved  to  be  robust. 

8.  RECOMMENDATIONS 

In  order  to  make  the  whole  design  even  simpler  more  returning  gains 
loops,  of  the  Simplified  State-Feedback  and  Estimator  autopilots,  must  be 
eliminated.  A  further  investigation  then  must  be  conducted  in  order  to 
examine  whether  the  performance  of  the  resulting  design  remain  unchange. 
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APPENDIX  A 

GEOMETRIC  AND  MASS  PROPERTIES  OF  MISSILE  CONFIGURATION 

Weight  (lbs) 

2525 

Maximum  Diameter  (in) 

24 

Length  (in) 

168 

1,0,  (Slug  Ft2) 

40 

Iyy  (Slug  Ft2) 

804 

Izz  (Slug  Ft2) 

810 

Center  of  Gravity  Distance  from  Nose  (in) 

100.8 

Reference  Length  (Ft) 

2 

Reference  Area  (Ft2) 

Tt 
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APPENDIX  B 


ANALOG  TO  DIGITAL   CONVERSION  OF  UNCOUPLED  PITCH  CHANNEL 


AUTOPILOT 
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APPENDIX  C 
PROGRAM  LOGIC  FOR  COMPUTING  CONTROL  LAW  K  VIA  ACKERMANN'S  FORMULA 

1 .  Read  In  A,  B,  T,  and  Nj,  the  number  of  states. 

2.  Comment:  first  we  will  read  In  the  desired  pole  location  in  the 
plane,convert  them  to  z-plane  polynomial  coefficients,  and  construct 
<x(A). 

3.  I  *  Identity  matrix,  NsxN$ 

4.  ALPHA   *   I 

5.  k  +    1 

6.  If  k>  Nj,  go  to  step  18. 

7.  Read  in  pole  location  k  as  a+jb. 

8.  If  b=0,  go  to  step  14. 

9.  A1  *  -2  exp(aT)cos(bT) 

10.  A2*  exp(2aT) 

11.  ALPHA    *    ALPHA  X  (A  X  A  +  AtA  +A2I) 

12.  k  <■    k  +  2 

13.  Go  to  step  6. 

14.  Ai   +    exp(aT) 

15.  ALPHA   <■    ALPHA  X  (A  -  A,  X  I) 

16.  k  *  k+  1 

17.  Go  to  step  6. 

18.  Comment:  now  we  construct  the  controllability  matrix. 

19.  C   +    I 
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20.  E    +    B 

21.  k  <+    1 

22.  If  k  >  N,,  go  to  step  28. 

23.  Comment:  replace  column  k  of  C  by  E. 

24.  CI  ;k]  +  E 

25.  k  +    k  ♦  I 

26.  E   <r  A  x  E 

27.  Go  to  step  22. 

28.  Comment:  now  solve  for  the  control  law,  first  form  eTn  as  the  last 

row  of  I 

29.  E  «■    UN,  j    1 

30.  Solve  DC  =  E  for  D. 

31.  K  =  DxALPHA 

32.  END 
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APPENDIX  D 


COMPUTER   PROGRAM   FOR  APPLICATION   OF   ACKERMANN'S  FORMULA 


FILE:      ACKERMAr,    «ATrIV  Al 

IJCfi 

INTEGER  NS,  K..IA,  IE,^,  I  ,J,  Nl 

REAL  PH!2C,2i,),GAl2C,2J).AI(2G,2O),ALPHA(20,20; 

REAL  A2        ,c  I  20, <:J)  .0  (  20,  2  0)  t  S120)  ,  oM<:  J, 2  J)  ,  oiU  23, 

REAL       A,o,T,     R,  ATEMV120,  20) 

REAL  m1  ,  TG!2J,2o)  ,  Til  23, 20)  ,  TSl  20,20)  ,^  120,20), 

C 

C  INr'UT    OF    NS=    NIUM3E.<    01-    STaTlS 

C  T=    SaMPL::     lift 

PRINT    ,        '     ENTER    NUMBER    OF    STATES.' 

READ,    NS 

PRINT,      •    ENTER    SAMPlE     TIME    • 

READ,     T 
C  INPUT    JF    PhI=A     MATkIX 

00    100     1=1, NS 

DC     101     J=l,.,S 

WRI  TE(  6  ,103  )      I,  J 
103       FORMAT  (  4X,«  ENTER    P  H(  •  ,  12  ,  •  ,  •  ,  12,  '  )  •  ) 

REAO.PH!  I, J  ) 
101       CONTINUE 
100       CONTINUE 

c        input   of   oa=3   mat:;  IX 

OQ    2  00       1  =  1  ,NS 

J  =  l 

WRI TE(6  ,105  )      I,  J 
105       FORMAT!  4X,'  enlTEn    GA(  •  ,  12  ,  '  ,  '  ,  12,  '  )  '  ) 

REAO.GA  (I ,J  ) 
200       CONTINUE 
C  BUILDING    THE    IDENTITY    MATRIX 

C  Ai_PHA=C=l 

DO     1    1  =  1,, NS 
DO    5    0  =  i,NS 
IF     (  I  .£  C.J)     GC    TO    300 
All  I  ,  J ) =o . 
ALPHA!  I  ,j)  =  0. 
C(  1  ,J)  =  0. 
GL    Tu    5 
300        All  I  ,  J  )  =1. 

ALPHA!  I  ,J)  =  i  . 
CI  I  .  J)  =  1- 
5  CONTINUE 

1  CONTINUE 

K  =  l 
16  IF     (K.GT.NS)     Gu    Tu     Id 

C 

C  INPUT    OF    OEjIrEC    PULEo    LOCATION 

C 

«RITE     (6,43) 

43  FORMAT!  •    ENTER    REAL    PART     OF     OtSIREO    PuLE    LOCATION.') 
RtAu,A 
HKITE!  6  ,44 ) 

44  FORMAT!  '     ENTER     IMAGINARY     PART    OF    OESI/.ED    PULE     LOCATION. 
READ.o 

C        BUILDING    ALPHA!  PHI) 

IF     «B     .  EJ.O)     GO     TO     14 

A  1= -2.* EXP!  «*T)  *CCS!d*T) 

A2=EXP ( 2.*A-T) 
C  EVALUATE     PH*PHl 

CALL    VMULFF      (  PH  ,PH  ,  NS  ,  NS  ,  NS  ,  2  0  ,20  ,  T  M  ,2  0,  I  CR  ) 


C 


DC    2    I  =  1,NS 

DU    o    J=1,nS 

TQ ( J , J)spH| I, j)*Al 
TU  I  ,  J)=AI  !  I,  J)*A2 
TS(  I  ,  J)  =  TM(  I,J)+TU!I,J)+TI  (I, J) 
6  CONTINUE 

2  CONTINUE 

CALL    VMULFF     ! ALPHA , T S , NS , NS , NS  ,20 , 2 0 ,A TEMP, 20 ,  1 1 R ) 
DO    600     1=1, NS 
DO    601     J=1,NS 
601        ALPHA!  I  ,j)  =  ATEr(P(I  ,  0  ) 
6C0       CONTINUE 


ACKOOO  10 

AC  r\OcO  2 0 
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AC*C042C 
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ACXO  j4dO 
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AC. <0C5 00 
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AC *00 5 40 
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ACK.00560 

ACO'J570 

ACKCU580 
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ACK0U600 
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ACN006VO 

ACK00700 

A; *007 10 

ACi\0O720 
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FILE:     ACKEkMAN    *ATFIw  Ai 


1<» 


701 
700 


13 

9 

12 

10 


28 
11 


C 
C 

SENTRY 


K=    iv+2 

GO    TC     16 

Al=EXPl  A-*T) 

00    3    I  =  1,NS 

00    7    J=1,NS 

T0(  I ,J)  =Ai  (  I  ,  J) *A1 

TI(  I ,J> =PH(  I  ,  JJ-TCI  I , J  J 

CCNT  INUE 

CALL    VMULFF  I  ALPhA,  T  I  ,  N3 ,  N  S  ,N  j  ,  20  ,  20  ,  A  7  EMP  ,  2  0 ,  I  E  ?. ) 

00    700     1  =  1,  .4S 

CO    701     J«L,,ilS 

ALPHA!  I  ,J)  =  ATEMrMI  ,J  J 

CCNT  INUE 

K=Kt-l 

GO    TO     16 

OG    9    1  =  1,NS 

J  =  l 

E(  1  ,  JJ  =  GA(  I  ,  J  ) 

K=l 

IF    (K.GT.N3)     GO    TG    28 

00    10    I=l,NS 

C(I  ,K)  =  EU.  1  ) 

CONTINUE 

K=K*1 

CALL    VMULFF      I  PH  ,  E.  N:, ,  NS  ,  1     . 2  J  ,  2 J , c, 20,  IEK i 

GO    TO    11 

DO    11    J=1,NS 

El     1  .  J  J  =AI  I  ,*StJ  i 

R*0. 

CALL    LG  INF     I  C  ,2  0  «.«S  » '>tS  »r\  ,  D  ,  2  0  ,3  »  *r\ ,  I  E;<  ) 


CALL    VMULFF 
CALL     VMULFF 


IE,J,1     ,NS,NS ,20 , 20, 3f ,20, IEK) 

IbF  ,^LPHA,  1     ,  NS.  f.S  ,2  0,20,  JK,  20,  IERJ 


PRINT,      »  CONTROL     GAIN       VECTQ* 

WRITE     (6,^5)        IUMi,JJ,J-liNSI 

FOR.HAT  (  oH  5  .4  ) 

RETURN 

END 
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AC i\0u9O 
aC*00950 
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ACK j1U60 
AC  kG  10  70 
AC. <0  10^)0 
AC/.01U90 
AC  is 01 100 
ACi\Gli  10 
AwK01120 
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APPENDIX  E 

ANALOG  TO   DIGITAL   CONVERSION   OF  COUPLED   PITCH  AND  ROLL   CHANNEL 
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APPENDIX  F 

ROBUSTNESS   DATA   FILE   FOR  THE   DISCRETE   STATE   FEEDBACK   DESIGNED 

AUTOPILOT 

FILE:     FEEDBACK    JATA  Al 
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FILE:     FEEDBACK     DATA 
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FILE:     FEED8AC*     OAT 
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APPENDIX  G 

ROBUSTNESS   DATA   FILE  FOR  THE   DISCRETE   ESTIMATOR  DESIGNED 

AUTOPILOT 
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